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I. OVERVIEW 



It has been nearly ten years since the ISGW model [1] was introduced [2,3] so it is not 
surprising that the heavy quark semileptonic landscape now looks very different. At that 
time, for both theoretical and experimental reasons, inclusive decays were the main focus of 
attention, and the ISGW model, which studied exclusive decays and approximated the inclu- 
sive semileptonic spectra by summing over resonant channels, was considered quite eccentric. 
Today, improvements in both theory and experiment have made exclusive semileptonic de- 
cays a main focus of attention. Such decays seem very likely to provide the most accurate 
determinations of the weak mixing angles V c b and V u b- They also provide excellent probes 
of hadronic structure via precision tests of Heavy Quark Symmetry (HQS) [4-7]. 

The ISGW model was in many respects a stepping-stone to Heavy Quark Symmetry: 
it is a model which respects the symmetry in the heavy quark limit near zero recoil. It 
also played a role in the discussion of the reliability of the free quark decay model (and its 
derivatives) for the endpoint region in b — > u semileptonic decay. Indeed, the model had its 
origin in that discussion, and was designed to provide the minimum reasonable prediction 
for the decay rate in this region for a fixed V u (,. In this paper we present an updated version 
of ISGW, which (with the permission of the ISGW authors) we call ISGW2 to emphasize 
that it is not a new model but rather an improved version of an old one [8]. The new features 
are described in detail in Section III, but briefly they are: 

1. Heavy Quark Symmetry constraints on the relations between form factors away from 
zero recoil are respected, 

2. Heavy Quark Symmetry constraints on the slopes of form factors near zero recoil are 
built in [9], 

3. the naive currents of the quark model are related to the full weak currents via the 
matching conditions of Heavy Quark Effective Theory (HQET) [6], 
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4. Heavy-Quark-Symmetry-breaking color magnetic interactions are included, whereas 
ISGW only included the symmetry-breaking due to the heavy quark kinetic energy, 

5. the ISGW prescription for connecting its quark model form factors to physical form 
factors is modified to be consistent with the constraints of Heavy Quark Symmetry 
breaking at order I/thq, 

6. relativistic corrections to the axial coupling constants (known to be important in the 
analogous coupling qa in neutron beta decay) are taken into account, and 

7. more realistic form factor shapes, based on the measured pion form factor, are em- 
ployed. 

The discovery of Heavy Quark Symmetry has not eliminated the need for models; it has 
rather provided a solid foundation for model-building and redefined the role that models 
should play. Consequently, an updated version of the ISGW model that incorporates the 
lessons of Heavy Quark Symmetry, and is designed with current usage in mind, seems very 
worthwhile. Among other roles, models should: 

1. provide predictions for the various universal form factors ("Isgur-Wise functions") of 
Heavy Quark Symmetry, 

2. provide predictions for the form factors governing b — > u, c — > s, c — > d, and s — > u 
transitions not directly governed by Heavy Quark Symmetry, and 

3. give estimates for the sizes of Heavy-Quark-Symmetry-breaking effects in the b — > c 
decays determining V c b, in the relations between b — ► u and c —>■ d matrix elements 
which can be used to determine V u b from exclusive semileptonic decays [4,10], and in 
the relation between c — » s and b — > s matrix elements which enter into the prediction 
of exclusive b — > decays [10]. 

In the next section we will give some of the background to ISGW and to the events 
leading up to ISGW2, as well as a quick review of the basic elements of the ISGW approach. 
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As already mentioned, Section III describes the new features of ISGW2 in detail. In Section 
IV we present our results. Section V discusses their implications for Heavy Quark Sym- 
metry, while Section VI compares our results to experiment. Section VII closes with a few 
comments. 
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II. BACKGROUND 
A. Some History 

In 1985, when the model that was eventually published as the ISGW model [1] was 
introduced [2,3], its intended use was very different from its present use. Moreover, much less 
was known about semileptonic b and c quark decays, both theoretically and experimentally. 
The ISGW2 model presented here is designed to update the earlier version to address both 
of these shortcomings. Ten years ago, the experimental study of the semileptonic decays of 
b and c quarks was in its infancy. In particular, for b quarks the main available data was 
on the inclusive lepton energy spectra for B — > X£i?£, generated by the quark level b — > clT>i 
and b — > utvt transitions. At that time the principal theoretical tool being used to analyze 
these spectra was the QCD-corrected parton model of ACCMM [11] and its relatives [12], 
with particular emphasis on extracting the Cabibbo-Kobayashi-Maskawa (CKM) [13] matrix 
elements V c b and V u t, from inclusive lepton spectra. Early fits to these spectra [14] near the 
b — > civn endpoint were leading to alarmingly small upper limits for the ratio \V u b/V c b\ 2 . 
Such results could of course simply be attributed to errors in the data. Alternatively, they 
could be taken as serious limits on V u b which would indicate a failure of the Standard Model 
scenario for CP violation. The ISGW model was introduced to explore a third possibility: 
that a partonic description of the b — > clvg and b — > uivi transitions in the endpoint region, 
where the lepton energy is near its maximum, might be deficient. The basic motivation for 
this concern arises from the observation that the highest energy leptons in these decays are 
asssociated with the production of the lowest-mass hadronic final states X in B — > X c £i>g 
and B — > X u ti>t respectively; the partonic description would only be expected to apply once 
the states X c and X u had masses above their respective resonance regions. 

We will discuss this issue in more detail below. We raise it at this point to recall that 
one of the main goals of the ISGW model was the production of an alternative description 
of the endpoint region which intentionally represented an extreme example of how little 
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b — > uii>i could show up in the endpoint region. The motivation was to illustrate the the- 
oretical uncertainty which should be reflected in upper limits on \V ub /V c i,f extracted from 
inclusive endpoint spectra and to thereby place more realistic constraints on Standard Model 
CP-violation scenarios. Along the path to this primary goal, the ISGW model produced a 
number of other results. In retrospect, the most important of these were probably concep- 
tual: much of the framework for Heavy Quark Symmetry [4-7] was presented in these early 
papers [1-3], including the vital role of the zero recoil point (where t = {pt + Pv e ) 2 is at its 
maximum value t m ), the insensitivity of B — > Dti>i and B — > D*£up transitions to mb/m c , 
and the role of D — > K£ + ue and D — > K*£ + i>e measurements in "tuning" exclusive models 
to be used for the extraction of V cb and V ub . ISGW also made a number of predictions. For 
example, ISGW was the first exclusive model to calculate rates to channels other than the 
pseudoscalar and vector ground states and consequently to predict that in both b — > c£p£ 
and c — > s£ + U£ decays the exclusive transitions B — > D, D* and D — > K, K* would dom- 
inate. This prediction (which is surprising since kinematically masses up to m B and m D , 
respectively, are allowed), now has a firm basis in theory [15,16, 4-7]. They also pointed 
out that in the nonrelativistic limit (applicable to such exotic processes as B c — > ipi&t), the 
weak transition form factors would be controlled by a set of universal functions given by 
the Fourier transforms of wave function overlaps and not by t-channel meson masses. This 
point has since been explored by many authors [17]. 

As mentioned in Section I, this update of ISGW has been prompted by a number of 
developments. The most fundamental of these is the discovery and development of Heavy 
Quark Symmetry [4-7]. In particular, the development of Heavy Quark Effective Theory 
[6] as a tool for systematically treating both the 1/rriQ and perturbative QCD corrections 
to the extreme Heavy Quark Symmetry limit has helped place models like ISGW in clear 
focus. HQET divides the calculation of current matrix elements into two steps: matching 
the currents of the full theory onto those of a low energy effective theory associated with 
some relatively light renormalization scale //, and then calculating matrix elements in the 
low energy effective theory. From this perspective, a quark model like ISGW or ISGW2 
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is presumed to be associated with a quark model scale /i qm ~ (9(1 GeV) where a valence 
constituent quark structure of hadrons dominates the physics. 

Since the constraints of Heavy Quark Symmetry for current matrix elements of the low 
energy effective theory are consequences of QCD, every model should display these results 
(including an allowed symmetry-breaking pattern) in the appropriate limit. In fact, in the 
low-recoil region where nonrelativistic dynamics apply, the ISGW model was already totally 
consistent with the Heavy Quark Symmetry limit. Adding the constraints of Heavy Quark 
Symmetry in ISGW2 nevertheless has significant impact. In high recoil b — > ctvn transitions, 
some ISGW form factors have missing functions of w = v ■ v' (v and v' are the four- velocities 
of the initial and final hadronic systems; this variable is called w after the origin of the 
name of this letter in, e.g., French) which are unity at zero recoil, e.g., the / form factor in 
B — > D*£ui is missing a factor of ^(1 + w) which goes to unity at w — 1. A related issue 
is embedded in the recoil dependence of the ISGW form factors. As discussed in ISGW, 
the slope of a quark model form factor consists of two terms: a normal "transition charge 
radius" term and a relativistic correction (of order l/m^mj in a Qi — > Qj current matrix 
element) which is outside of the scope of a nonrelativistic quark model. ISGW posited that 
such relativistic effects could be taken into account in an approximate way by replacing all 
factors of (t m — t) appearing in their nonrelativistic formulas for form factors by K~ 2 {t m — i), 
where k is the ratio of the nonrelativistic charge radius to the true charge radius. Heavy 
Quark Symmetry [9] tells us that this prescription (while fortuitously close numerically in 
the cases to which it was applied) is incorrect; the symmetry moreover dictates the correct 
result in the heavy quark limit. This result, to be described below, is adopted in ISGW2. 
Consideration of the allowed pattern of HQS-breaking at order 1/mg also has an impact. 
Among other effects, it requires a change in the ISGW prescription for relating the form 
factors of the weak binding limit calculated here to physical form factors. Although such 
modifications to ISGW are only strictly required near the heavy quark limit, ISGW2 adopts 
the usual constituent quark model stance of treating all constituent quarks like heavy quarks, 
so the same changes are made, e.g., to c — > s transitions. 
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QCD also demands that the matrix elements of a low energy effective theory like the 
quark model be corrected by the matching conditions which map them onto the matrix 
elements of the full theory. At the level of the currents of the two theories, these matching 
conditions take the generic form 

J£ = C,,J£ + f AJ£ + —^—SjJji + J$ . (1) 

" L Qj " L Qi 

In ISGW2, we explicitly calculate the l/mQ j and l/mQ i corrections in the quark model, so 
only the factor Cjj mapping the naive vector (Qj7 M Qj) and axial vector (QjT^sQi) cur- 
rents of the quark model onto the true currents [Qj^^Qi and QjT^TsQ?) and the expansion 
in terms of the new naive currents appearing in AJ^ in order a s /ir are needed. We will give 
these matching factors in Section III. A below. 

There are other reasons why an update of the ISGW model is warranted. In the pe- 
riod since the publication of ISGW, its role in providing a very conservative upper limit on 
\Vub/V c b\ 2 has become antiquated; ISGW2 attempts to modernize ISGW so that its predic- 
tions become best estimates rather than most conservative estimates. Consider, for example, 
the curve in Figure 1 showing the ISGW form factor F W (Q 2 ) with Gaussian wavefunctions. 
The charge radius of the pion was used to determine the value of the parameter k which 
in turn determines the rate of decrease of F n (Q 2 ) shown. Thus, instead of choosing a value 
that provided a best global fit to the data over the whole kinematic range applicable to the 
b — > uivg transition, ISGW chose a value that fits at low Q 2 but, as a result of its unrealistic 
Gaussian form, underestimates F V {Q 2 ) at high Q 2 . This choice was driven by the ISGW 
goal of providing a minimum rate for B — > X u ti>t in the endpoint region. In ISGW2 we 
attempt a more realistic description of the recoil dependence of all form factors. 

There have also been important experimental developments since 1985! In B decays 
[19], the inclusive spectra near the endpoint region show a definite B — > X u iv l excess [20], 
although, for the reasons already mentioned, the resulting value of V u b is unclear. The decays 
B — > D£u£ and B — > D*£i?f have both been measured [21] in sufficient detail to extract the 
CKM matrix element V c b with some confidence since the observed features of these decays 
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are consistent with the expectations of Heavy Quark Symmetry. Preliminary evidence for 
B — > D**£u e decays (here D** represents non-D or D* decays) has been reported and searches 
have begun for the exclusive b — > ulvi processes B — > pti>n and B — > uive [22] . In D decays 
[23], where V sc is known, D — ► K£ + u e and £> — > K*£ + is e decays have been measured [24] in 
sufficient detail to extract the four c — > s form factors contributing in the limit — > 0, and 
rather tight limits onD -> K**£ + i>e have been set. In all cases the experimental results are 
qualitatively consistent with the predictions of ISGW (despite some initial indications to the 
contrary [25]); indeed, all results to date are consistent with ISGW within its anticipated 
"quark model accuracy" of predicting matrix elements to ±25%. However, in the spirit of 
"tuning" the quark model to higher accuracy, in ISGW2 we have taken note of a substantial 
failure of ISGW to predict the magnitude of the S-wave axial vector form factor / in D — > 
K*£ + Vi decay. In the quark model, this form factor is analogous to g& in neutron beta 
decay, where experiment is about 25% below the quark model prediction of 5/3; the data on 
/ indicate that it is also smaller than the quark model prediction. There is a very natural 
explanation for the g& discrepancy within the quark model [26,27]: the matrix elements 
of the space components of the axial current in a relativistic S-wave spinor are reduced in 
proportion to the probability of lower components in that spinor. We accordingly build this 
relativistic correction factor into ISGW2. 

With the predictions of Heavy Quark Symmetry to facilitate the extraction of V cb and 
V u b from the data, one of the main uses of models has shifted from predicting form factors 
to predicting the deviations of form factors, or relations between form factors, from the 
predictions of Heavy Quark Symmetry. In view of this changing role, we implement one 
further elaboration of ISGW in ISGW2: we consider the effects of hyperfine interactions 
on meson wavefunctions. ISGW already naturally took into account the other 1/rriQ effect 
in HQET [6], the heavy quark kinetic energy, so this addition to the model completes the 
parallel with the most general symmetry-breaking effects allowed. As we will see, the u qa 
effect" and these hyperfine interactions, in concert with matching corrections, eliminate the 
problem with the D — > K*£ + vt form factor /. 
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To summarize: ISGW2 is an updated version of ISGW designed to make "best estimates" 
within the context of a constituent quark model that fully respects Heavy Quark Symmetry. 

B. A Review of the Foundations of the ISGW Model 

In Section III we will describe in detail the new features which we incorporate in ISGW2. 
Here, we review the basic ideas and methods of the ISGW model. 

ISGW breaks the problem of computing a current matrix element of a transition from 
a state H of mass, momentum, and spin m,p,s to H' with m',p',s' into kinematical and 
dynamical parts. It first makes the usual mechanical Lorentz-invariant decomposition of 
the matrix element into Lorentz tensors and invariant form factors f\ (i = 1,2, ...N) which 
depend only on the four momentum transfer variable (t m — t) where t — (p' — p) 2 and where 
t m = {m! — m) 2 is the maximum momentum transfer. The variable (t m — t) is used since it is 
zero at the "zero recoil point" where H' is left at rest in the rest frame of H; the importance 
of momentum transfers near t m will be made clear below. 

It should be noted that any specification of the functions fi(t m — t) leads to a Lorentz 
invariant description of these weak decay processes. In this sense ISGW is not a nonrela- 
tivistic approximation. It is, however, a nonrelativistic estimate of the intercepts fi(0) and 
"charge radii" r,- L = [6 ^^ ]^ (or more generally the shapes) of the Lorentz invariant form 
factors fi(t m — t). These estimates are made by noting that there is a one-to-one correspon- 
dence between the fi and a partial wave expansion of the (if' Ij^O^if ) matrix elements. For 
example, if H and H' are pseudoscalars P and P', then 

(p\ P ')\v»(o)\p(p)) = fr(p+p7+f- p (p-pr ■ (2) 

This decay also has two partial wave amplitudes. In the rest frame of P, V°(0)\P(0)) is still 
a pseudoscalar, and so creates P' in an S-wave; V(Q)\P(0)) is an axial vector and so must 
create P' in a P-wave. Thus m(/f p + f?' p ) + E'(f p ' p - f p ' p ) and p'(f p ' p - f p ' p ) are 
proportional to the rest frame S-wave and P-wave amplitudes, respectively. 
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A vital element of the ISGW model is that each partial wave amplitude is calculable in 
the nonrelativistic limit; the one-to-one correspondence with the fi then allows a calculation 
of each in this limit (i.e., none of the fi are intrinsically relativistic in character). The ISGW 
model therefore calculates the fi in a limit in which the model would in principle be exact, 
and then extrapolates these exact formulas to the physical regime. It should be noted that 
the nonrelativistic limit requires more than p/m and p'/m' being small. It also requires 
that the internal motion of the constituents of H and H' be nonrelativistic. One of the 
essential assumptions of the ISGW model is that such "mock meson" form factors fi which 
are derived in the approximation that m u , ma, and m s are large compared to Aqcd can be 
extrapolated down to their actual constituent quark masses to estimate the fi. 

In the heavy quark world in which the ISGW formulas would be exact in principle, the 
low-lying mesons would all be simple quarkonia. ISGW is therefore necessarily a model for 
matrix elements between resonances, i.e., it does not directly address the issue of whether 
semileptonic meson decays are resonance dominated. The original ISGW paper argues that 
nonresonant contributions are likely to be small (their absence is correlated with the known 
success of the narrow resonance approximation), and there is some evidence from the data 
for this prediction. Nevertheless, the issue remains a hotly debated one. Note that this 
debate is relevant to our updating of the resonant matrix elements only once we use them 
to estimate the inclusive rates, e.g., those in the B — > X u ti>t endpoint region. At that point 
we will discuss this issue in more depth. 

The semileptonic decays of the be meson B c via the C7 M (1 — 75)6 current provide a 
good illustration of a system in which the ISGW model would in principle be an excellent 
approximation. Both B c and the low-lying states of the cc system can be reasonably well- 
described as nonrelativistic, and matrix elements like (ip(p's')\A u \B c (p)) can be accurately 
calculated in the frame where p = for small p' as atomic-physics- type wavefunction overlap 
integrals. This is the essence of the ISGW method. However, serious model dependence can 
occur when these matrix elements are extrapolated to large recoils; moreover, it occurs even 
at small recoil when any quark mass is extrapolated down to the constituent masses of the 
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u, d, or s quarks. Amongst the issues which must be addressed when a light quark plays a 
role are: 

1. the quarkonium approximation: It is a fundamental tenet of the constituent quark 
model that, up to "small" corrections which arise from pair creation leading to reso- 
nance widths, systems containing a light quark can still be treated as quarkonia. I.e., 
extra qq pairs and the gluonic degrees of freedom do not have to be introduced ex- 
plicitly. ISGW adopts this approximation. In addition, for simplicity it uses harmonic 
oscillator wavefunctions to approximate the true quarkonium wavef unctions. 

2. the weak binding approximation: For heavy quarkonia, the quark masses and energies 
are approximately equal and as a result the hadron mass is approximately the sum 
of the constituent quark masses. Once a quark becomes light, the failure of this 
approximation and other associated complexities make even a low velocity boost of 
the wavefunction problematic. Moreover, when both constituents are light, there can 
be a great discrepancy between the mass of a hadron and the sum of the masses of its 
constituents; this in turn leads to ambiguities in the extrapolation of the nonrelativistic 
formulas. ISGW adopts a specific prescription for dealing with such ambiguities. 

3. relativistic corrections: Even if the extrapolation of the nonrelativistic ISGW formulas 
were straightforward, they would still suffer from their failure to incorporate impor- 
tant relativistic physics. A simple example is the charge radius r\: in general such a 
radius will receive both nonrelativistic contributions with a scale controlled by the ra- 
dius of the quarkonium wavefunctions and relativistic contributions controlled by the 
Compton wavelengths of the participating quarks. (The latter contributions are them- 
selves of several types: relativistic corrections from the quarkonium wave equation, 
field-theory-induced pair creation effects, etc. Such effects are simply lost in the non- 
relativistic limit since r ~ 1/p ^> 1/m.) Of course the form factor intercepts will also 
receive relativistic corrections: generically, we can write fi(0) = /™ r (0)[l + 0(p/m)]. 
ISGW invokes the empirical success of the nonrelativistic quark model in assuming 
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that /™ r (0) will be a reasonable approximation to /j(0). Finally, there is of course no 
guarantee that fi(t m — t) won't contain (t m — t) dependence that is intrinsically rela- 
tivists. For example, an additional "kinematic" factor of [1 + *™ ~* ] would be "seen" 
as unity by a nonrelativistic calculation. 

All of these shortcomings, and others left unmentioned, make it surprising that the 
nonrelativistic constituent quark model works as well as it does. It may be that its successes 
are based on one crucial fact: that "it is better to have the right degrees of freedom moving 
at the wrong speed than the wrong degrees of freedom moving at the right speed" [28]. Given 
that the quark model would be correct if all the quarks were heavy quarks, its utility may 
reside in its ability to parameterize the evolution of the properties of these correct degrees 
of freedom from heavy to light systems. 
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III. ISGW2: THE NEW FEATURES 



As already repeatedly emphasized, ISGW2 is not a new model: it is a slightly improved 
version of ISGW [1]. In this section we describe one by one the differences between ISGW 
and ISGW2. For closely related studies of the marriage of the quark model with the physics 
of Heavy Quark Symmetry, see the work cited in Ref . [7] . 



A. the constraints of Heavy Quark Symmetry 

Although ISGW is completely consistent with the constraints of Heavy Quark Symmetry 
at maximum recoil t m (or w — 1) in the symmetry limit, HQS also determines various 
aspects of the behavior of the form factors at finite recoil and at nonleading order in the 
1/rriQ expansion. For example, the six form factors of B — > Divg and B — > D*£i?£ are 
required by Heavy Quark Symmetry to have, in the low energy effective theory, the form 
[4-7] 

/+ + /- = /+-/-= 9 = f = a- - 6+ = £(w) (3) 



a + + a_ = (4) 

where £(w) is the Isgur-Wise function (we have adopted conventions for defining the HQS 
form factors which lead to these simple forms; see Section V for explicit formulas relating 
these to the usual ISGW form factors). In the heavy quark limit, the ISGW model 
respects all of these constraints at all w except for that on /: it gives / = (j^)£(u>), 
corresponding to the nonrelativistic approximation 1+w ~ 2. Such effects, which correspond 
to v 2 /c 2 corrections to a leading nonrelativistic prediction, lie outside of the dynamical 
framework of ISGW, but are easily appended to the model (see, e.g., the second of Refs. [7]). 
Using eq. (3) and the corresponding results of Ref. [29] on the B — > D\lv^ B — > D^ 2 £u£, 
B — > D^ 2 £p£, and B — > D^iv^ decays, the required modifications to be incorporated into 
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ISGW2 are easily enumerated. They are all listed explicitly in Appendix C. The simplest 
example is the axial vector form factor / of eq. (B15) of ISGW, which, as the form factor 
corresponding to / in eq. (3), picks up an additional factor of \{1 + w) in the HQS limit 
(additional nonleading effects in the 1/mg expansion will be described below). In addition 
to these modifications, Heavy Quark Symmetry tells us that in heavy quark systems the 
eigenstates with J p = 1 + are not the L — S coupled states 3 Pi and 1 P±, but rather the j — j 
coupled states P^ 2 and P^ 2 with s\ l = | + and | + , respectively [30,29]. We therefore also 
list the new form factors appropriate to semileptonic decays to such excited P-wave mesons 
in Appendix C. 

In addition to these constraints of Heavy Quark Symmetry on the matrix elements of 
the low energy effective theory, HQET prescribes how to match these matrix elements onto 
matrix elements of the full theory, as already mentioned above. The matching of a generic 
form factor fj^ of type a associated with the underlying Qi — > Qj transition can be written 
in the form [31] 



7T 



(5) 



where the are unity for /+ + /_, /+ — /_, g, f, and a + — a_, and zero for a + + a_. Here 
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Nf the number of active flavors below the scale to; (four for i — b) and Nf the number below 
rrij (three for j = c). In contrast, the radiative correction functions j3j"\w) multiplying a s /ir 
(evaluated at a scale /iji intermediate between m; and rrij which we take to be the geometric 
mean ix 3 i = (mjirii)^) are a— dependent. 

The j3j"\w) associated with each of six form factors /+, /_, g, f, a + — a_, and a + + a_ 
are known. At w — 1 they are simply 



where 



and 



with 



^ ++/ - ) (l)=7ii-^-* 



(10) 
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We deviate from the use of these matching conditions only in case of transitions between light 
(u, d, and s) quarks. Since, as described in Appendix A, we assume that a s "freezes out" 
at the quark model scale /i qrn , the "renormalization group improved" matching conditions, 
as embodied in the Cji factor, are inappropriate for such transitions. For them we use the 
expansion of Cji to lowest order in a s , i.e. we resort to "lowest order matching". 

In principle the ftffi are functions of w, but this w— dependence is predicted [31] to be so 
weak relative to uncertainties in the w-dependence associated with nonperturbative effects 
that we ignore it here. (This dependence would, for example, correspond to a change in the 
predicted rates for the exclusive b — > c decays of the order 1% if it could be distinguished 
from the w— dependence in the preasymptotic nonperturbative Isgur-Wise functions.) On 
the other hand, there is nontrivial w— dependence contained in the factor 
w near 1, 



a s (f-lqm ) 



: for 



a L (w) 



cl- 2 - 8 ,)ln 
3 \33-2N' f J 



a s (fi 



qm J 



a s (mj) 



(w - 1) (19) 



We will make use of this factor below. 

The preceeding HQS-and HQET-induced modifications to ISGW are consequences which 
emerge from considerations of the heavy quark limit. There are additional modifications 
which arise from restrictions on the form of I/ttlq corrections to this limit. From the most 
general form of l/m b and l/m c corrections to the B — > Dlv f and B — > D*£u e form factors 
(see Section V), it is possible to resolve an ambiguity in the procedure for relating the 
form factors of the ISGW weak binding nonrelativistic calculation to physical from factors. 
Such a calculation in principle only determines form factors up to factors like (mj//mii-) n 
where m# is a physical hadron mass and rhu is the sum of its constituent quarks' masses. 
HQET resolves this ambiguity in a pleasing way: it specifies that the form factors ff m 
being calculated in such a quark model are (up to order 1/mg) the dimensionless form 
factors fi of the heavy quark limit which expand matrix elements in Lorentz invariants 
using the heavy quark four velocities t> M , t>' M and not the f\ which expand them in terms 
of their momenta. (As far as we can determine, it is purely by accident that the ISGW 
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notation for the weak binding form factors coincides with the notation for the HQS form 
factors). Such considerations in addition demand that the conventional form factors /; be 
obtained from the ff m by mass scaling factors which differ from the physical masses by at 
most 1/rriQ effects. In ISGW2 we choose to resolve this residual ambiguity by using the 
hyperfine-averaged physical masses rfiH 1 H 2 = ^fZi = r ^ l H 2 °f a HQS spin doublet of hadrons 
H 1 and H 2 to relate the ff m to the /j. (For the s,z — \ ground state doublet this mass is 

just rfi V p = \m v + \m P - in general m Sl = (j^) m j= Sl +\ + (2^+1) m j=s t -\ for an H Q S 
multiplet with light degrees of freedom having spin si). As in ISGW, we use the /; to 
compute all rates; this may be viewed as a residual model-dependent choice of certain l/trig 
terms, and illustrates very clearly how HQS and HQET have reduced the model-dependence 
of the results of ISGW2 relative to ISGW. 

There are two additional but clearly related elements to the correspondence between the 
ff m and the The ff m are functions of w which is the weak-binding variable analogous to 
the physical variable w. In passing to the physical form factors which depend on t m — t 
we identify 

yj _ 1 = tm ~ t (20) 

2mp Q m Xq 

for a transition Pq — > X q ivg induced by an underlying Q — > q£i/i transition. In addition, we 
correct all form factors for relativistic terms proportional to t m — t that are required by the 
form of 1/wiq corrections: see Section V and Appendix C for details. 

We next note that consideration of sum rules for the Qi — > Qj£&e transition in the heavy 
quark limit provides a constraint on the slope of the Isgur-Wise function £(w). If we define 
p 2 by the expansion 

£M = l-p>-l) + --- (21) 
about w — 1, then [9] in the heavy quark limit 

P 2 = \ + Pl d + Herf (22) 
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As will be seen below, the structure-dependent term p 2 sd dominates for a weakly bound 
system. The Ap 2 ert term is that appearing in eq. (19) from the w-dependence of the 
matching factors. The \ represents a relativistic correction to the nonrelativistic limit. It 
corresponds to the relativistic correction found in Ref. [27]; its generalization to systems 
with different spins is discussed in Ref. [32]. The results of Ref. [27] actually dictate the 
subleading (order 1/mg and l/m q ) corrections to the | of the heavy quark limit. In terms 
of a conventional charge radius r 2 defined by 

/(f) = f(t m )[l- l -r 2 (t rn -t) + ---} (23) 

in an expansion of the generic Pq — > X q form factor / around t — t m , the relation corre- 
sponding to eq. (22) with subleading terms from Ref. [27] included is 

Cts \fJ>qm) 

_ a s (m q ) _ 

where, for a ground state harmonic oscillator wavefunction (see Appendix A), 



y.2 — i- t 2 H - ( — I In 

Am Q m q wf m PQ m Xq \ 33- 2N' f J 



(24) 



2 g,p 

T ' w -f = 2m PQ m Xq Pl x ■ (25) 

The terms of eq. (24) are associated in order with the terms of eq. (22). Indeed, r 2 j 
in eq. (25) is the transition matrix element of the square of the interquark separation 
between Pq and X q ; it would be four times the squared charge radius of the pion in the 
case where Pq = X q = ir. Since m sp /j3p X ~ msp/Vwf ^ 1 in the nonrelativistic limit, the 
| is indeed a "relativistic correction", as stated earlier. However, in the constituent quark 
model m sp /p w f ~ 1, so it could be a very significant "correction"! ISGW recognized the 
generic possibility of 1/mgm, corrections to r 2 and accordingly introduced a "relativistic 
correction factor" k to compensate for them: they took f(t m — t) — > /(^i^), corresponding 
to enlarging r by a factor kT 1 . It is now clear that the 4m ^ m term in eq. (24) is a well- 
defined and necessary relativistic kinematic correction which should be added to the r 2 j 
term of a nonrelativistic model. In ISGW2 we note that this required correction is actually 
sufficient to achieve the same empirical effect as the multiplicative factor k of ISGW, which 
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was fit to the low t pion charge form factor. While relativistic dynamics missing from the 
constituent quark model might in principle still affect r^, we assume that such effects were 
on the whole subsumed into the quark model once its parameters were chosen to give a 
good description of the meson spectra. Thus, in ISGW2 we drop the ad hoc k factor in 
favor of the use of eq. (24). This has the additional bonus of making ISGW2 consistent 
with the dynamical constraints of the Bjorken sum rule [9]: in the heavy quark limit, the 
"structure-dependent terms" in eq. (22) are determined by the amplitudes to excite final 
states X q with sj e = | + and | + . Finally, we note that the App ert term vanishes for decays 
to s, d, and u quarks since their masses are already below the quark model scale /i qm where 
the running coupling constant has been assumed to saturate (see Appendix A). As a result, 
it only comes into play for b — > c transitions. 

In addition to this improvement in the way we deal with the slopes of the form factors, in 
ISGW2 we also abandon the gaussian form factors of ISGW, which are unrealistic at large 
recoils. This modification is described in Section III.C below. 

B. some relativistic corrections to the quark model 

As stressed in Section II, the ISGW model was introduced to illuminate some basic issues 
surrounding semileptonic decays. It therefore used the simplest possible version of the quark 
model capable of addressing these issues. It is, however, known that the predictive accuracy 
of the naive nonrelativistic quark model can be substantially improved by considering various 
relativistic corrections to that model. One of the simplest such corrections occurs in the 
matrix elements of the axial current. The naive nonrelativistic quark model predicts that 
qa — | in neutron beta decay. However, it has been known for twenty years that when the 
constituent quarks are given realistic momenta, g& is reduced by a factor of 1 — ^Pi ower , 
where Pi OW er is the probability of lower components in the quark spinors [26]. By taking this 
effect into account, most models [26,27] obtain values of g A about 25% smaller than |, close 
to its observed value of 1.257 ±0.003 [33]. For the Pq — > V q axial vector S-wave form factor 
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(called / in ISGW) the correction factor is 



C ' = J v.*J d ^M-2Er\ t 1 - 3(E q + rn q )(E Q + m J l-2E^\ 

(26) 



where the 0's are S-wave momentum space wavef unctions, Ei = (p 2 + m 2 )^ , and Jy q p Q = 
J d 3 p(f>y q (p)(f)p Q (p) . For a heavy quark transition in the heavy quark limit Cf — 1, but for 
a light quark transition, the analog of eq. (26) for n — > p would give roughly the required 
reduction of qa- In ISGW2 we adopt this correction factor as being at the least a reasonable 
interpolation between these two extremes. The correction factors resulting from eq.(26) 
using the masses and wavefunctions of Appendix A are given in Table I. 

Two potential deficiencies of this approach should be noted. There is in the first place 
no reason to suppose that there are not other more dynamical effects which renormalize the 
matrix elements of the light axial quark currents: the effect taken into account by eq. (26) 
should be only part of the story [34] . In addition, it is not clear that only the S-wave form 
factor / will be affected by relativistic corrections. We nevertheless take this as the simplest 
working hypothesis, and assume that the effective constituent quark mass subsumes other 
relativistic corrections as it does for quark model magnetic moments [27]. 

A second class of relativistic corrections to the quark model appears in the wavefunc- 
tions themselves. For simplicity, ISGW ignored the effect of relativistic corrections to the 
effective interquark potential. In particular, although quark model hyperfine interactions 
are responsible for the B* — B, D* — D, K* — K, and p — n splittings, their effects on the 
wavefunctions were not taken into account. (In HQET [6], this origin of the B* — B and 
D* — D splittings can be given a firm foundation via the o lxv G lxv /2rriQ operator appearing 
at order 1/toq in the heavy quark expansion. The quark model assumes the continuing 
relevance of this mechanism for light quarks as well). For our purposes, the net effect is that 
pseudoscalar and vector particles of a given flavor are no longer characterized by the same 
wavefunction parameter /3s (see Table II of ISGW). An update of this Table which takes 
into account this splitting is given in Appendix A. Given that both hyperfine and spin-orbit 
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effects in P-wave mesons are empirically very weak, we ignore such effects. 



C. more realistic form factors 



As mentioned above, ISGW used the gaussian form factors generated by their highly 
truncated harmonic oscillator basis; moreover, they used them out to relativistic recoils. 
Here we attempt a more accurate parameterization of the form factors which will have a 
more realistic behaviour at large (t m — t) by making the replacement 



cxp 



K2 
6' 



rl f (tm-t) - l + -^r 2 (t m -t) (27) 



-N 



where r 2 is given by eq. (24). In eq. (27), N = 2 + n + n' where n and n' are the harmonic 
oscillator quantum numbers of the initial and final wavefunctions (i.e., N = 2 for S-wave 
to S-wave, N = 3 for S-wave to P-wave, N = 4 for S-wave to S'-wave, etc.). These form 
factors all have the charge radii dictated by the quark model in the nonrelativistic limit, 
approach the gaussian form factors of the harmonic oscillator model as N — > 00, but provide 
a much better global fit to the pion form factor (see Fig. 1). In fact, with eq. (24) we predict 
< r 2 . >^= 0.61 fm, in satisfactory agreement with the observed [18] value of 0.71 ± 0.02 
fm. Since the Q 2 range covered by this figure corresponds to a (t rn — t) range that covers 
the recoils available in the semileptonic decays we treat here, we adopt the substitution of 
eq. (27) for all our decays. We emphasize that these substitutions should be viewed as low 
energy parameterizations of the form factors and not as appropriate descriptions of their 
analytic or high (t m — t) forms. 
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IV. RESULTS AND DISCUSSION 



The formulas we require to predict semileptonic form factors and rates may all be ob- 
tained from ISGW [1] (supplemented by formulas given in refs. [25] and [29] for decays in 
which the lepton mass is not negligible) by making the few simple modifications described 
in the text. The required changes are described explicitly in Appendix C. To calculate rates 
we insert into these formulas the constituent quark masses and j3- values from Tables Al and 
A2 of Appendix A. 

We now present our results, organized by the underlying quark decay and arranged in 
order of increasing spectator quark mass. We will compare these results to the predictions 
of Heavy Quark Symmetry in Section V and to experiment in Section VI. 

A. b^clui 

These decays are generally the most stable predictions of our model, and those that are 
underwritten by Heavy Quark Symmetry are the most reliable. All states contain a heavy 
quark and the available recoil is limited, reducing the sensitivity to form factor slopes. Since 
the b and c quarks are not only heavy but also have a modest mass difference, the Shifman- 
Voloshin limit [15,16] is also relevant to the decays with a light spectator and thus provides 
a simple explanation for why the electron spectral shape is very similar to that of the free 
quark decay model despite dominance by the ground state pseudoscalar and vector final 
states. 

Our results for B— >X C( ^eu e are shown in Figure 2; the partial widths are given in Table II. 
This decay is dominated by the pseudoscalar and vector meson final states, which contribute 
29% and 61% of the total semileptonic rate respectively. Our absolute prediction for the 
inclusive decay rate for B^>X c gtv e is T = 4.06 x 10 13 | V c b\ 2 sec -1 , about the same as the 
ISGW result. The approximate validity of the SV limit gives an electron spectral shape 
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very similar to the free quark model despite the dominance by the ground states of X c . Our 
predicted form factors are also close to those of the heavy quark limit; a detailed discussion 
of this limit will be given in Section V. 

2. B s ^X cS £u e 

The small difference in the m d and m s constituent quark masses on the m b or m c scale 
results in B s — >X cS ei> e decays behaving in a very similar fashion to the previous case. This 
may be seen in both Figure 3 and Table II where our results are displayed. As expected, 
there is a small increase (~ 5%) in the total fraction of the IP and 2S states compared to 
B decay since the SV limit holds here to a slightly reduced degree. Our absolute prediction 
for the inclusive decay B s ^X cs eu e is Y = 3.90 x 1 13 1 V c b\ 2 sec -1 , slightly smaller than the 
previous case. 

3. B c ^>X C c£&e 

This decay is different from the preceding two cases in several ways. As the spectator 
quark is no longer light, both the parent and daughter mesons are approximately non- 
relativistic and are thus appropriately described by our model. The results, which should 
therefore be quite reliable, are shown in Figure 4 and Table II. They are still reminiscent 
of the previous results with lighter spectators even though the spectator approximation 
prediction that the inclusive semileptonic decay rates should be equal fails by about 25% in 
going from B d to B s to B c deacys. The contributions from the pseudoscalar and vector final 
states are, however, reduced (to 29% and 47%, respectively) as expected from the spectator 
arguments given in the ISGW papers and from the inapplicability of the SV limit. 

The measurement of the slopes of the form factors for these decays would provide an 
interesting test of the arguments made in Refs. [1,2,17] that naive dispersion relations for 
these slopes will fail. These systems are predicted to have charge radii determined by their 
Bohr radii ~ [ f^ l" 1 while dispersion relations would lead one to believe (unless one 
were very careful [17]) that the charge radii will be of order (mj + m c ) -1 . Discussing the 
possibility of studying these states may not be completely far-fetched: there are suggestions 
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[37,38] for experiments to observe them. 



B. c^sl + ui 

The decays induced by the quark level process c^s£ + vi are dominated by the ground 
state pseudoscalar and vector daughter mesons. This is a consequence of the low available 
recoil momentum which has little probability of producing excited states. These decays are 
nevertheless not expected to be as accurately described as the b^c case since the s quark 
is too light for Heavy Quark Symmetry to apply. 

Our results for D^X sQ e + v e are shown in Figure 5; the partial widths are given in 
Table III. This decay is predicted to be almost totally dominated by the pseudoscalar and 
vector meson final states, which contribute 63% and 34% of the total semileptonic rate 
respectively. Our absolute prediction for the rate of the inclusive decay D^X sU e + u e is 
T = 0.17 x 10 12 |Vs C | 2 sec -1 , down about 10% from ISGW. This decrease arises from an 
increase of the K rate of 18% and a decrease in the K* rate of 41% which dramatically alter 
the ISGW K* j K ratio; for details, see Section VI below. Our total predicted width is about 
one half that of a simple free quark model using our constituent quark masses. 

It is amusing to note that while neither the heavy quark nor SV limits should be appli- 
cable here, they both seem to have strong residual influences on this decay. The comparison 
of our predicted form factors with those of Heavy Quark Symmetry assuming that s is a 
heavy quark will be given below. 

2.D s ^X s - s i+vi 

While m u and m s are very similar on the scale of the charm quark mass, they are 
noticeably different on the scale of the daughter quark mass m s . As a result we do not 
expect as strong a similarity between D and D s decays as that which existed between B and 
B s decays. In addition, non-ideal mixing in the pseudoscalar sector of D s decays leads to 
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a very different spectral shape due to the comparatively low mass of the 77. This low mass 
gives a much higher electron endpoint than the corresponding free quark decay endpoint, 
as may be seen in Figure 6. The fractions of the semileptonic rate going to 77, 77', and are 
31%, 26%, and 40%, respectively, with the distribution of rate to the 77 and rj' sensitive to 
the assumed pseudoscalar mixing angle of —20°, but with the sum relatively insensitive (see 
Table III). The IP and 2S contributions are once again predicted to be small: only 3% of 
the semileptonic width. The ratio of rates for D^K*e + u e and D s ^xf)e + is e is 

T{D^K*e+u e ) ~ °- 84 ' (28) 

down about 20% from the ISGW value [39] of 1.02. This decrease is mainly due to hyperfine 
interaction effects and the new rh prescription of Section III. A. We note before leaving these 
decays that the D s inclusive semileptonic decay rate is itself down by more than 25% from 
D decay. This substantial failure of the spectator approximation will be discussed in Section 
VI. 

Z.B c -+Xrf+v t 

Our results for B c -^X s ie + u e are shown in Figure 7. The explicit partial widths are 
given in Table III. Not surprisingly, the extreme mass of the spectator in this case results 
in a spectrum that is very different from the naive free quark spectrum. It is dominated 
by decays to the pseudoscalar (43%) and vector (55%), as the available energy is small. 
Recoil effects are very small due to the large daughter mass. The softening of the lepton 
spectra expected due to the high spectator mass is pronounced as is the reduction of the 
inclusive rate. Our absolute rate T(B C — >X s ie + u e ) = 0.50 x 10 11 |V^ C | 2 sec -1 is less than a 
third that of D— >X s{i e + z/ e , corresponding to an even more dramatic failure of the spectator 
approximation. 

The ratio of B c decay via b^c to c^s decay is 



T(B c ^X sl e+v e ) 
T(B c ^X c5 e + u e ) 



0.0014 



Vsc 



2 



(29) 
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for \Vcb\ ~ 0.04 and \V SC \ ~ I. Thus, amusingly, B c semileptonic decays will be roughly 
evenly split between the two very different quark level processes b^c and c^s. 

C. c^dl+vt 

We now consider the Cabibbo-suppressed decays involving the quark level process 
c^dt^vn which are also predicted to be dominated by the ground state pseudoscalar and 
vector final states. These decays have taken on a new importance since the realization that 
their measured form factors can be related to the form factors of b— m via Heavy Quark Sym- 
metry. As indicated in Section V below, these relations should eventually lead to accurate 
model-independent determinations of |Kt&|. 

\.D Q ^X m i + v t and D + ^X dd i + v, 

Our results for D°^X du e + is e and D + ^X dd e + u e are shown in Figures 8a) and b) re- 
spectively. The partial widths are given in Tables IV and V, respectively. D — >X du e + v e is 
dominated by n and p final states which contribute 63% and 31% of the total respectively, 
compared to 43% and 52% in ISGW. This shift in relative probability comes mainly from 
a substantial decrease in the p rate. However, there is also a sizeable 5% rate predicted to 
the J p = 1 + P-wave states. The longitudinal to transverse ratio for the p is 0.67. The D + 
decays look somewhat different as the final states now include both the 1 = and 1=1 
neutral states. Note that T(D + — >X dd e + is e )/T(D ^X d ue + v e ) = 0.92, which is mostly due to 
the effects of the 77 and rj' channels. We also note in passing that Cabibbo-forbidden decays 
are predicted to represent approximately 5% of D° decays and 4% of D + decays. 

2.D s + ^X ds £ + u e 

D s ^X ds e + v e decays of Figure 9 and Table IV are again dominated by the pseudoscalar 
and vector ground states which contribute 60% and 29% of the total resonant semileptonic 
rate with a 10% contribution from the J p = 1 + P-wave states. The absolute rate is almost a 
factor of three times smaller than the free quark decay rate, and the lepton spectrum is much 
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softer. Note that these Cabibbo-forbidden decays are predicted to contribute approximately 
6% of the inclusive D s decay rate. 

D. b^uivt 

We now consider the decays corresponding to the quark level process b^-uiue. These 
decays are very important in the determination of the \V u b\ matrix element, which is itself 
important for CP violation in the Standard Model. Here large recoils are available; as 
a result we are not surprised to find large contributions from the IP and 2S states in 
our lepton spectra. In ISGW the decays to the radially excited pseudoscalars n 1 ^ were 
explicitly checked to confirm that the calculation would converge to the inclusive rate in the 
appropriate limit. As in ISGW, however, ISGW2 only sums over the low-lying resonances 
and so for b — > u decays it can be used as a model for the inclusive spectrum only in the 
endpoint region. This point is discussed at greater length in Section VI. 

1. B° -► X ud £v e and B~ X uU iu t 

We consider both the decays B '^X u jez/ e and B~ ^X uU eu e which are shown in Figure 
10a) and b), respectively. Detailed partial widths are given in Table VI. As in ISGW, there 
are large contributions from the IP and 2S states. On comparing with the results of ISGW, 
one sees that our more realistic form factors have increased the rate to the sum of the rates 
to the IS*, IP, and 2S states by about 25% and transferred some of the rate from the heavier 
to the lighter states. ISGW2 therefore predicts a somewhat hardened endpoint spectrum 
relative to ISGW. The change in individual exclusive rates is most pronounced for the pion, 
which has increased by about a factor of four over the ISGW result. As discussed above, 
and as is apparent from Fig. 1, ISGW was designed to produce a conservative estimate of 
the endpoint rate. The effect on the n rate is uncharacteristic since it vanishes for kinematic 
reasons at zero recoil where the ISGW form factor is nearly equal to ours, and grows into 
the high recoil region where their form factor is far below ours and the measured pion form 
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factor. The consequent large uncertainty in the tt rate is compounded by the potential effect 
(to be discussed below) of the nearby B* pole. In contrast, the ISGW2 rate to the p is only 
about 70% larger than that of ISGW. A similar increase is obtained for the total rate to the 
IP states, while the rate to the radial excitations of the ir and the p decreased by almost a 
factor of three. We also note that the ratio T L /T T for B^pev e has remained equal to the 
ISGW value of 0.30, even though this value is sensitive to the method used to treat large 
recoils and the axial current form factors. 

As in ISGW, we have included all resonance states with masses < 1.7 GeV, which implies 
that our lepton spectra are complete for lepton energies greater than about 2.4 GeV. In this 
region our spectra are considerably softer than that of the free quark decay. Since our sum 
over exclusive channels is incomplete, we cannot quote a total rate, although the treatment 
of the pion radial excitations described in ISGW suggests that it will be within a factor of 
two of our free quark rate of 1.28 x 10 14 |Kb| 2 sec -1 . 

2. B s ^X u - s tv t 

Our results for B s — >X uS ev e are shown in Figure 11, with the explicit partial widths in 
Table VII. As expected, this decay is very similar to that of the 5-meson. There is, however, 
a noticeable softening of the spectrum due to the heavier spectator. 

3. B c — ► X u5 lvt 

Our results for B c — *X U cev e are shown in Figure 12. The explicit partial widths are given 
in Table VII. This decay is similar to the other b^uev e decays. However, the softening of 
the spectrum due to increased spectator mass is much more pronounced, as is the shifting 
of probability to states with masses above those of our calculation. 
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V. COMPARISON TO HEAVY QUARK SYMMETRY 

In this section we compare our results to those of Heavy Quark Symmetry, which pro- 
vides model-independent predictions for some aspects of the weak hadronic matrix elements 
presented here. While these model-independent predictions are very interesting theoreti- 
cally, the sizes of the corrections to this limit may restrict its validity to a limited number 
of processes, or to a small region of phase space. One way to estimate the effects of such 
Aqcd/^<? corrections to these limiting predictions is to compute these corrections in a model 
such as ISGW2. ISWG2 is in fact most reliable precisely at the key zero-recoil point of Heavy 
Quark Symmetry, and indeed our form factors reduce to those required when taken to the 
symmetry limit. Away from this limit our results constitute model-dependent predictions 
for the effects of the finite quark masses. Estimates of such corrections from other hadronic 
models and from quenched lattice QCD have also been made [40]. 

The predictions of Heavy Quark Symmetry for the decays P Q ^P q £p£ and P Q ^V q £p£ were 
first worked out in Refs. [4] . The relationship between our form factors and those of Heavy 
Quark Symmetry described by eqs. (3) and (4) which are defined in terms of four-velocity 
variables is 

2 V V m P* V mp Q / 2 \\ m P q V mp Q / 

2 v V m p q V mp Q / 2 v V mp * V mp Q / 

g = 2 v /m PQ m Vq g 

f ~_ f(l + w)- 1 
^m PQ m Vq 

m 2 P 

a + + a_ = (a + + a_) 

v /m PQ m Vq 

a + - a_ = - v /m PQ m Vq (a+ - a_) . 

Recall from Section III. A that since these relations involve the physical hadron masses m H , 
and not the hyperfine-averaged masses rh H , these /, differ from the ff m calculated directly 
at the quark model level by l/rrig terms. 
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(30) 

(31) 
(32) 
(33) 

(34) 
(35) 



The l/m q and 1/mg corrections to these predictions have been considered by various 
authors [41-45] . In particular, the first of Refs. [45] gives a general form for such corrections 
which is, as we will see, particularly suited to our quark model, namely 

/+ -1 + eM (36) 



/- J_(A 



+ P4H (37) 



1 + tA" + — P2(w) + — (piH - p 4 H) (38) 



f(tu) 2/x+ m g m Q 



f , A (w — 1\ 1 . . 1 / . . w — 1. 

1 + — — +— p 2 (w) + ( Pl (w) —Pa(w)\ (39) 



£(u>) 2/i + \w + lj m q m Q \ w + 1 

fa, + a) 



- (^TT) ^ (K — (w+ l)p 3 (w) + p,{w)) (40) 



£(w) (w + 1) 

(a+ - a_) , A/«)-ll 1 \ 1 , 

• -1 + - ——— + —+— UH-/03H r-rP4(^, 



£(w) 2 l^to + lm, m Q/ ™ ? V w + 1 

+ — ( Pl (w) - p 4 (w)) (41) 



where A is a constant and l/p± = l/m q + 1/mg. The inclusion of these effects thus results 
in the appearance of four additional unknown functions p n {w) with unknown normalizations 
(although it can be shown that pi(l) = ^2(1) = 0). Alternative parameterizations have also 
been given; we will comment on one of these below. 

We can map onto our results by expanding them to leading order in \jm q and l/m Q ; we 
find 

« m -i-M!0 (42) 



t P+ 
f qm m sp 
£ = 2/x_ 



(43) 



£T = 1 + !^££ _ _ (44) 

/ qm 1 i? y (w) i? P (w) , /w-l\m sp 



= 1 _ ^^^^ _ 11£1 ^ + „, . ^ 

4 m g mg Vw + l/2/i + 



~qm . ~qm 



m g (l + u>) 



(46) 



and 
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g^ m - a qm = 1 + msp _ m sp _ R v (w) _ R P (w) 
£ 2fi + m q (l + w) m q m Q 



Here 



*w<«0 i + ^-o — (48) 

with 5 sp = mpQspVQsp — m Q as mg — > oo and with the A's parameterizing the approach of 
the /9's to the heavy quark limit via 

with Ap sp = /c sp — |/i sp and \v sp = k sp + \h sp . Here k sp describes the perturbation of (3q sp 
in a heavy quark meson with heavy quark Q and a spectator sp due to the heavy quark 
kinetic energy (which is spin-independent) and h sp is the analogous perturbation due to the 
residual hyperfme interaction of Q and sp. From Table A2 one can see that kd ^ +0.14 
GeV, h d ~ +0.36 GeV, k s ~ +0.26 GeV, and h s ~ +0.50 GeV. Using the measured masses 
and the constituent quark masses of Table Al, and correcting for the residual heavy quark 
kinetic energy, one can estimate that 5d ~ 0.09 GeV, and 5 S ~ 0.17 GeV. 
This decomposition allows us to identify 

A = m sp (50) 

Pl (w) = -R P (w) (51) 

p 2 (w) = -R v (w) (52) 

PsH=0 (53) 

p 4 H=0 , (54) 

from which one can easily see that the predicted corrections to the Heavy Quark Symmetry 
limit are all of modest size. This assessment is made quantitative by Table VIII for b^c and 
c^s decays with an up or down spectator. In terms of the common ipi, ifj 2 , ^3, ip+ param- 
eterization [44] of 1/m.Q and l/m q corrections, these results are ip\ = —(Rp + 3i?y/2A)£, 
ip2 — 0, ^3 — ~{Rv — -Rp/4A)£, and ip + = —(w — l/w + 1)£. Note that, as required, ipi 
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respects heavy quark spin symmetry, while ips is responsible for breaking it. Conclusions 
similar to ours (couched in terms of the ^-parameterization) have been reached previously 
in the quark model [7]. 

As an aside, let us note that if we focus on B — > D*£u e transitions alone, and assume 
only that p^{w) = p±{w) = 0, we may define a "preasymptotic Isgur-Wise function" Cd*b{ w ) 
such that 
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1 + 
1 + 



A 



f 

£d*b(w) 
(a+ + a_) 
£d*b(w) 

(a+ -a_) = l+ A 



A /iy - 1 



2/x + Vu> + 1 
A 



m q {w + r 



A 



(55) 
(56) 
(57) 
(58) 

£ D * B (w) ^ ' 2/i+ m 9 (w + 1) 
Under this assumption, therefore, the predictions of Heavy Quark Symmetry to order l/m,Q 
and l/m q can be described by one unknown parameter A known to be approximately m sp 
and the unknown shape of the "preasymptotic Isgur-Wise function" (,d*b( w ) which retains 
its normalization to unity at w — 1. We also note that the measured slope of this function 
is predicted by our model to be 



P 2 d*b - 0.74 



(59) 



which value includes a contribution 



A 2 16 , 

A Ppert ~^l ln 



qm J 



0.13 . 



(60) 



a s {m c ) 

(The numerical value of App ert is quite uncertain: it depends on the leading logarithmic 
expansion in m c / ' p qm and on the assumption that a s (p qm ), where p qm is the "quark model 
scale", is the "frozen-out" value a s = 0.6 from Appendix A. We accordingly assign a theo- 
retical error of ±0.05 to it.) 

Another important set of predictions [10] of Heavy Quark Symmetry are those which 
relate, e.g., the form factors of B — > pti>n to those of D — > p£ + ve- These predictions could 
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play a vital role in determining V u b if corrections to the symmetry limit are not too severe. 
With ISGW2 we can check these relations. For example, in the ideal symmetry limit one 
should have 

« S K)-^ K) ^=^ = «.kr w ^^ (6i) 

2^/m p m B 2^/m p m D 
where p p and vp are the four momentum of the p and the four velocity of the decaying 
meson P and ai(m,Q) is given by eq. (7) with Nf appropriate to rriQ. (Note that m p has 
no special significance in these formulas: we are simply using it to create dimensionless 
quantities. Also note that we have removed the known quark mass dependence of the 
leading logarithmic matching condition, but not attempted to remove the mass dependence 
contained in the a s /n corrections since, while relatively weak given that a s (fi u b) — Q: s (^dc), 
it is model-dependent.) We find, e.g., that at zero recoil 

a s fa)" a;K) , 1 P) = 0.49 (62) 
2 v /m p m s 

aJm,)- a '^ f ^ P{mp) = 0.45. (63) 

One also expects 

2a s (m b )- a ^^fT-mg B ^ = 2a s {m c )- a ^ jW~^g D ^ (64) 
while our model predicts (once again at zero recoil) 

2a s (m b y a ^^-mg B - p (m p ) = 1.16 (65) 

2a s {m c )~ ai{mc) ^/m p m D g D ^ p (m p ) = 1.15. (66) 

The form factors /+ and a + are more complex since it is the combinations /+ ± /_ and 
a + ± a_ which obey simple scaling relations. However, since the objects which scale are 
yfij%(f+ + /-)» \f^(f+ - /-)> m B^f^( a + + a_), and ^m B m p (a + - a-), in the heavy 
quark limit /_ = — /+ and a_ = — a + so that in fact the simple scaling laws 
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y rriB y mo 

and 

2a s (m b ) - ai {mb) ^m p m B a 1 l^ p = 2a s (m c ) ~ ai (mc) y/m p m D a+~" p (68) 
emerge. Our model in fact gives 



a s (m b )- a ' (m »\ ^f^(m n ) = 0.68 (69) 

V m B 



a.W^g/r'K) = 0.66 (70) 



and 



2a s (m b ) ~ a/ (m&) y /m p m B a+^ p (m p ) = -0.66 (71) 

2a s (m c )- a/(mc) v / m p m ^ a +" P ( m p) = -°- 60 - ( 72 ) 

We conclude that our model strongly supports the conclusion that l/rnq effects will not 
obscure the extraction of V ub for exclusive B decays via the scaling relations of Heavy Quark 
Symmetry so that the proposal [10] to do so appears to be sound. (It should be noted 
that in the case of f+^ n and /+^ 7r , our quark model contributions at zero recoil must be 
supplemented by the B* and D* pole terms [46] , respectively, before they may be compared 
to experiment. These pole terms carry with them large but known I/tuq effects related to 
the smallness of relative to the B* — B and D* — D hyperfine splittings.) 

Similar conclusions follow for the validity of relations between c — > s and b — > s matrix 
elements which enter into the prediction of exclusive b — > 57 decays. 
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VI. COMPARISON TO EXPERIMENT 
A. Magnetic Dipole Decays 

Magnetic dipole decays of mesons like ui — > 7r7, K* — > Kj, and ip ~^ Vcl proceed 
through a transition magnetic dipole moment form factor which is precisely analogous to 
the vector current form factor g in weak decays of ground state pseudoscalar mesons to 
ground state vector mesons. The ability of our model to describe such decays is therefore 
relevant to the reliability of the model for the weak decays which are the focus of this 
paper. For the transition magnetic dipole moment [ipy = hvp underlying the transition 
V — > P7 (or P — > V7 when it is energetically allowed), theory (experiment [33]) gives, in 
units of the nucleon magneton, fi np = 0.52 (0.69 ± 0.04), fj, nu) = 1.56 (2.19 ± 0.09), = 
0.07 (0.13±0.01),/i w = 2.16 (1.77±0.17),/v = 0.68 (0.57±0.07), ^ = 0.61 (0.66±0.02), 
Aw = 1.53 (1.20 ± 0.08), AW = 0.58 (0.42 ± 0.04), ^ = -0.94 (|^| < 1.8), h k +k*+ = 
0.95 (0.79 ± 0.03), fi K o K *o = -1.27 (-0.98 ± 0.26), and AVv = °- 76 (°- 55 ± °- 12 )- As 
in the main calculations we have taken the pseudoscalar mixing angle here to be —20°; we 
have also assumed that the vector mixing angle is 39°. 

We conclude from this comparison that the quark model will probably be able to predict 
the form factor g with the typical quark model accuracy of ±25% for transitions involving 
light quarks. Since Heavy Quark Symmetry guarantees that our formulas for g will be 
correct in the heavy quark limit, this should be an upper bound to the probable error in 
such predictions. 

B. K^-kIvh 

Although the form factors for these decays are usually referred to the 577(3) symmetry 
normalization point t = 0, we prefer to refer them to the point t = t m where Heavy Quark 
Symmetry will develop. We find that f+(t m ) = 1.04 and /_//+ = —0.28. The latter is in 
reasonable agreement with the measured value [33] though there is a substantial uncertainty 
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since K ± decay gives /_//+ = — 0.35±0.15 while K° L decay gives — 0.11±0.09. Our equation 
for /+ is consistent with the Ademollo-Gatto theorem [47] which protects /+ from substantial 
deviations from unity. Our prediction for f+(t) can be compared with the "standard" 
[48,49] used to extract V us from these decays. If we convert to the linearized form f+(t) = 
/ + (0)[1 + \r^ K t\, then we predict /+(0) = 0.93 and r\ K = 0.48 fm versus the "standard" 
[49] /+(0) = 0.97 ±0.01 and r\ K = 0.53 fm corresponding to K* pole dominance. The best 
current fit value to this transition radius gives [33] r\ K = 0.59 ± 0.02 fm. 



C. Meson Decays through b^clvi 

Our results for semileptonic meson decays involving the quark level decay b^clvg were 
given in Section IV. A. Their relatively low recoil and heavy quark masses provide a theo- 
retical stability that makes them our most reliable predictions. 

Reviews of the experimental status of semileptonic B meson decays can be found in Refs. 
[19]. From the measured rate (here we use the latest CLEO result [21]) 

T(B^D*£u e ) = 2.99 ± 0.39 x 10 10 sec _1 , (73) 

and our prediction T(B^D*£u e ) = 2.48 x 10 13 |V r c6 | 2 sec~ 1 we obtain 

\V cb \ = 0.035 ±0.002. (74) 

The measured rate for B^Divi is 

Y(B^D£v l ) = 1.3 ± 0.3 x 10 10 sec" 1 . (75) 

Using our predicted rate of Y{B^Dlv^) = 1.19 x 10 13 | V cb \ 2 sec" 1 implies that 

| V cb \ = 0.033 ±0.004. (76) 

The consistency between eqs. (74) and (76) of course means that the model correctly predicts 
the ratio of the rates to D and D*. However, these determinations of \V cb \ depend on the 
prediction of the recoil dependence of the relevant form factors and so have a theoretical 
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error of order 10%. A comparison with data near zero recoil using Heavy Quark Symmetry, 
as has become standard, remains the reliable way to determine \V c b\. 

The measurements of Yl/Yt for B~^D*eu e are quite sensitive to the relative importance 
of the f,g, and a + form factors. Experiment gives 



0.85 ±0.45 (ARGUS [21]) 

0.83 ± 0.33 ± 0.13 (CLEO [21]) (77) 
1.24 ±0.16 (CLEO [50]) 

consistent with our prediction of 1.04 (versus 0.97 for ISGW). Furthermore, the predictions 
of both ISGW and ISGW2 for the q 2 dependence of B°^D* + £u e agrees reasonably well 
with the measured results: see Refs. [19]. In particular, we predict that the slope of the 
preasymptotic Isgur-Wise function £d*b( w ) wm be p 2 D ,^ = 0.74 ± 0.05 (see the text below 
eq. (60) for an explanation of the theoretical error), while the latest fits to the data (see the 
last of Refs. [21]) give for the closely related quantity p 2 the value 0.84 ± 0.14. (The ISGW 
prediction was 0.69.) Table IX shows our predictions for the individual form factors in terms 
of the HQS form factors defined in Section V. It also compares them with the predictions 
of ISGW, Heavy Quark Symmetry, and HQET (with matching but no 1/mg corrections). 
This comparison illustrates the relatively model-independent nature of these predictions. A 
recent measurement [50] gives for B — > D*£u£ decay 

j = 0.031 ± 0.009(stat) ± 0.004(syst) GeV~ 2 (78) 

°± = -0.015 ± O.OOQ(stat) ± 0.003(syst) GeV~ 2 (79) 

in reasonable agreement with our predictions of 0.030 GeV~ 2 and —0.024 GeV~ 2 , respec- 
tively. 

Both CLEO and ARGUS currently find indications that the D and D* final states ac- 
count for much less than all of the semileptonic decay width of the B meson. We predict 
that these final states account for ~ 90% of the total rate to the states included in our 
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calculation. If confirmed, these experimental results may indicate that non-resonant pro- 
cesses are important, or, perhaps, that we have underestimated the effects of the IP and 2S 
states. Note that the Bjorken sum rule [9] requires that the rate missing from the D and 
D* channels be approximately proportional to p 2 . Thus doubling the missing rate would 
require doubling p 2 , in apparent contradiction to the existing agreement between theory and 
experiment described above. 

Finally we note that the decays of the Bd, B s , and B c sequence show a marked departure 
from the spectator approximation in which their inclusive semileptonic decay rates would 
all be equal. This phenomenon, which is more pronounced in the c — > s decays, is addressed 
in the next subsection. 

D. Meson Decays through c^s£ + V£ 

The quark level decays c^s£ + ui are at this time better measured than the b^c£i>£ 
decays. They also provide a greater challenge for our model since in these decays Heavy 
Quark Symmetry does not guarantee the success of the leading approximation to their 
form factors: strange quarks do not qualify as heavy quarks! Note that since the CKM 
matrix element \V SC \ may be related to \V u d\ via the unitarity of the CKM matrix, direct 
measurements of the form factors can be made. The experimental status of weak charmed 
meson decays was recently reviewed in Refs. [23]. 

Averaging over measurements [24] and using isospin gives [23] 

T(D — > Ktvf) = 9.0 ± 0.5 x 10 10 sec -1 (80) 

which compares reasonably well to our prediction of 

T(D^K£ + u e ) = 10.0 x 10 10 sec" 1 . (81) 

In this decay one can also measure the pole mass for the /+ form factor assuming a monopole 
shape. CLEO obtains M^ lc = 2.00±0.12±0.18 GeV, consistent with earlier experiments but 
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with smaller errors. This mass corresponds to a transition radius r KD = —^2- = 0.24 ± 0.03 

^pole 

fm compared to our prediction of 0.22 fm. The data cannot currently distinguish between 
the common choices (monopole, dipole, exponential) for the shape of this form factor as 
the available range of t m — t is limited and all these shapes give an approximately linear 
dependence over this range. 

Assuming the measured form factor, the rate may be transformed [23] into a measurement 

of 

f+(t m ) = 1.42 ±0.25 (82) 

(or equivalently /+(0) = 0.75 ± 0.03). We predict f+(t m ) = 1.23 (or equivalently /+(0) = 
0.85 using our predicted t dependence and 0.80 using our form factor with the central 
experimental value of r KD . 

As an aside, we would like to explain why such form factor measurements should be 
referred to t — t m and not t — 0. Heavy Quark Symmetry establishes that heavy to 
light transition form factors are all related in the region of t m [4,10], i.e., are independent 
of the heavy quark mass mg as mg — > oo when scaled by an appropriate power of tuq. 
Measurements near t m are therefore determinations of universal transition form factors (up 
to 1/rriQ corrections). Form factors at t — are, in contrast, "random numbers" since they 
are the product of the universal transition amplitudes relevant at t m and a complicated 
dynamical function which depends on the microscopic details of the high momentum tails of 
the initial and final state wavef unctions. This is because t — corresponds to a final state 

9 9 

X recoiling with maximum three momentum \p x \ = 2 m P — m ^ ne res ^ fr ame °f Pq- This 
momentum increases with uiq so that t — form factors are ever- decreasing functions of 
m Q- 

The D meson semileptonic decay to the K* final state has been the subject of much 
interest. An early measurement found a value for Yl/Yt approximately two times larger 
than expected while the ratio of vector to pseudoscalar branching ratios was about one half 
what was expected from many models. Attempts were made [25] to accomodate these results 
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within the ISGW model by allowing for the theoretical uncertainties inherent to the quark 
model (±20%). It was found that the model could accommodate the vector to pseudoscalar 
ratio but not the IY/IV ratio with such variations. The current experimental situation is 
more precise with at least two independent measurements of each quantity. In addition to 
the above quantities, measurements of the form factors themselves have now been made. 

The averaged experimental measurements of Mark III, CLEO, E691, ARGUS, E653, and 
WA82 are [23] 

r(D^K*e + u e ) , . 

rWW) - M7± "' (83) 

T(D~^K*e+u e ) = 5.1 ± 0.5 x 10 10 sec" 1 (84) 

and 

^ = 1.15 ±0.17 . (85) 
It 

This compares reasonably well with our model values of 

^^'=0-54 (86) 

T(D^K*e + u e ) = 5.4 x 10 10 sec" 1 (87) 

and 

^ = 0.94 . (88) 

As anticipated in Ref. [25], agreement with the data relative to ISGW has come about 
via a modest shift in the form factor /. In fact, four different effects contribute: the matching 
conditions lower f(t m ) by 11%, Cf from Table I lowers it by about 7%, the wavefunction 
mismatch induced by hyperfine effects (see Table A2) decreases f(t m ) by another 7%, while 
the new factor of |(1 ± w) raises the average of the amplitude over the Dalitz plot by about 
7%. For this decay the form factors themselves have been determined. The comparison 
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of our predictions to the measured results [23] are given in Table X. Before leaving these 
decays, we note that (if we subtract our predicted Cabibbo-suppressed rate) the inclusive 
Cabibbo-allowed D semileptonic decay rate is measured [33] to be (16.2 ± 1.5) x 10 10 sec~ 1 . 
This compares favorably with our prediction of 15.8 x 10 10 sec _1 . 
Our predictions for the analogous form factors for D s ^(f)e + v e are 

f(t m ) = +2.03 GeV 

g(t m ) = +0.52 GeV- 1 (89) 
a+(t m ) = -0.29 GeV- 1 . 
These results are in reasonably good agreement with recent measurements [51] which give 
(see Ref. [50]) g(t m )/f(t m ) = (+0.20 ± 0.07) GeV' 2 and a + (t m )/f(t m ) = (-0.21 ± 
0.05) GeV~ 2 to be compared with our predictions for these ratios of +0.26 GeV~ 2 and 
-0.14 GeV- 2 , respectively. CLEO [51] also quotes T(D S -> r)e + v e )/T(D s -> 0e+z/ e ) = 
1.7 ± 0.4 and F(D S -> v( 'e + v e ) /T(D S -> 4>e + v e ) = 0.7 ± 0.2 to be compared with our predic- 
tions of 0.8 and 0.7, respectively, for a pseudoscalar mixing angle of —20°, and 1.2 and 0.5, 
respectively, for —10°. 

Ref. [51] also presents an extraction of the rate T(D S — > 4>e + v e ) based, among other 
things, on the assumption that the inclusive semileptonic decay rates of the D and D s are 
equal. This assumption would appear to be justified on the basis of recent work on the I/tuq 
expansion of inclusive heavy quark semileptonic decays [52-55]. However, it is inconsistent 
with the results quoted here, which predict that Y(D S — > Xe + v e ) is 27% smaller than T(D — > 
Xe + u e ), largely as a consequence of the restricted phase space in the r/' decay of the D s . We 
speculate that the unexpectedly [55] large corrections we predict arise from an inapplicability 
of the assumptions under which the strong version of the results of Ref. [52] were derived: 
since these decays (along with those induced by the b — > c transition) are dominated by the 
lowest few resonances, the spectral decomposition of the decay is imperfectly described by 
the smooth partonic spectral function. (As explained by the authors of Ref. [52], this is 
analogous to R in e + e~ annihilation being smooth and well-approximated by its partonic 
value only well above a threshold. We note that heavy quark semileptonic decays can be 
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deceptive in regard to when they are "well above a threshold" because, while the recoil 
mass can kinematically run up to the mass of the decaying quark, the hadronic spectrum in 
Q — > qtvi in fact cuts off at a recoil mass-squared only of order kQCDij^Q—m^) 1 /ttiq « ttiq 
above threshold in the decay of a Qd meson.) With the assumptions made, Ref. [51] obtains 
T(D S — > 4>e + v e ) = (4.4 ± 0.7) x 10 10 sec" 1 ; if their assumptions are modified to correspond 
to our predictions from Table III for the ratio of inclusive rates and for the degree to which 
T(D S — > (rj + rf + 0)e + z/ e ) saturates T(D S — >■ Xe + u e ), this extracted rate would be changed 
to (3.5 ± 0.5) x 10 10 sec" 1 . These results are both roughly consistent with our prediction 
that T{D S -> 0e+z/ e ) = 4.6 x 10 10 sec" 1 . 



E. Meson Decays through c^d£ + U£ 



The Cabibbo-suppressed charmed meson decays via the quark- level process c — > d£ + ui 
have taken on an enhanced importance recently. As described in Section V, Heavy Quark 
Symmetry relates the form factors for such decays near t — t m to their analogues induced by 
the crucial b — > uivg processes. In the short term the better measured c — > s£ + vg processes, 
combined with SU (3) flavor symmetry, can substitute for these decays, but precision deter- 
minations of |T4b| w iH probably require accurate determinations of the Cabibbo-suppressed 
form factors. 

Experimental studies of such decays have begun. The decays _D°^7r~e + z/ e and D + — > 
7r°e + z/ e have been measured by Mark III and CLEO II, respectively [24]. Using the value of 
iKd/VscI ~ 0.227 ± 0.003 which follows from CKM unitarity [33], their quoted results can 
be translated into the form 

J + ~ = 1-17 ±0.19 (90) 

where a pole model for the t-dependence of the form factors has been assumed in the deter- 
mination of the numerical factor, and we have averaged the results of the two experiments. 
Our model predicts the value 0.71 for this ratio. 
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F. Meson Decays through b^uivi 



Our results for semileptonic meson decays involving the quark level decay process b^utvn 
were presented above. Only the decays with a light spectator have been observed. Both 
CLEO [20] and ARGUS [20] have observed leptons in the 2.4-2.6 GeV energy range that 
can be populated only by leptons from a b^ulvi process. In addition, searches for exclusive 
modes like B — > irive, B — > ojtv^ and B — > pti>i have begun. 

The extraction of |14&| from these data is based on kinematics. For B mesons produced 
at the T(AS) resonance, decays via the quark level process b^clvt have a maximum lepton 
energy of 2.4 GeV/c, while the leptons from a b—>ulve process may have energies up to 2.6 
GeV/c. Consequently, these inclusive decay processes can be unravelled in the endpoint 
region of the lepton spectrum. 

As previously mentioned, the physics of this endpoint region has been the subject of 
intense discussion [56]. The ISGW papers met with strong criticism by many who argued 
that its treatment of the endpoint region was inconsistent with the parton model. This issue 
has recently been clarified in favor of ISGW by rigorous 1/mg expansions of the inclusive 
rate. The zeroth order argument was given in Ref. [57]. It is shown there how, in a b — > uive 
transition, the zeroth order lepton spectrum is controlled by quark level kinematics. The 
key observation is that decays to low mass hadronic final states (in this approximation, the 
hadronic mass is just the invariant mass of the recoiling u quark and the noninteracting 
spectator quark) only populate the high t (low recoil) region of their Dalitz plot which 
therefore cuts off their electron spectrum at the quark level endpoint energy. Thus while 
from kinematics alone such decays might have produced electrons with energies out to the 
physical (i.e. hadronic) endpoint, they do not for dynamical reasons. Conversely, high mass 
hadronic final states produce electrons out to their kinematic endpoint, but the highest 
mass hadrons have an endpoint which exactly coincides with the quark level endpoint in 
zeroth order. Recent work [52,53] has demonstrated that this picture is the beginning of a 
rigorous expansion of inclusive decays, and that the 1/mg corrections have exactly 
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the character anticipated by ISGW and Ref. [57] and continued in this work. 

It has very recently been speculated that, within the 1/mg expansion, even the endpoint 
region is amenable to treatment via an operator product analysis [54]. If, as indicated by Fig. 
10(a), this region is really dominated by a few resonances (mainly the p, a±, and 61), then 
this analysis may not apply. Thus while for b — > u decays, in contrast to c — > s and b — > c 
decays, the hadronic spectral function over most of the Dalitz plot will be well-approximated 
by the partonic spectral function (corresponding to strong applicability of the results of Ref. 
[52]), the endpoint region is only dual to the partonic spectral function in an average sense. 

The experimental analysis of this data requires simultaneous fits to the b—*ctvi and 
b^uivi inclusive spectra combined with the measured continuum backgrounds. Conse- 
quently, it is not possible to simply convert the values determined by various experiments 
which have used the ISGW model to the modified version of the model presented here, since 
the results are dependent upon the shape as well as the integrated inclusive rate. Since our 
spectrum is considerably harder than that of ISGW, it seems clear that the ISGW value 
of |Kib/K&| w iH decrease when reanalyzed. However, the change seems unlikely to be very 
large. While the rate to the p, which is most important in the extreme endpoint, has in- 
creased by 70%, the total rate to the states we consider has only increased by 23%. Given 
that these rates are proportional to |Vub| 2 , the decrease in (K^/V^l itself seems likely to be 
less than 25%. 

We would like to caution against interpreting this decrease, which brings ISGW into 
better agreement with other models of the endpoint region [58], as leading to a more reliable 
value for |T4&| from the inclusive spectrum. In the first place, it is a mistake to use models 
for this region which consider only the n and p final states: the endpoint region is clearly 
going to be populated by many more states than these. This exclusion leaves only ACCMM 
[11] and ISGW as potentially realistic models for this region. However, we would continue 
to stress, in spite of the real improvements of ISGW2, that our theoretical errors here are of 
order ±50%. The recent clarification [53-55] of the status of the ACCMM calculation [11] 
in this region suggests that it should be assigned a very substantial theoretical error as well. 
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VII. CONCLUSIONS 



We would argue that ISGW was already a good model for heavy meson semileptonic 
decay and that with the improvements added here ISGW2 is an even better model for 
this sector. ISGW2 behaves correctly in the Heavy Quark Symmetry and Shifman-Voloshin 
limits, including lowest order corrections to these limits. In taking into account the leading 
corrections to the Heavy Quark Symmetry limit, ISGW2 adds physics to ISGW which 
corresponds to that demanded by Heavy Quark Effective Theory. These corrections are 
implemented with a well-known, and well-tuned, model of quark dynamics. 

In order to extend the range of validity of the model and to include all relevant physics, 
ISGW2 also adds to ISGW other effects. These extensions have improved our agreement 
with the experimental data. For example, the mesonic decay rate for P—^VtV^ where P and 
V are pseudoscalar and vector mesons, respectively, is sensitive to the S-wave axial current 
form factor. This form factor probably receives sizable relativistic corrections (of order 10%) 
which we have attempted to take into account. The motivation for such extensions comes 
not only from first principles: in this case, such a correction is needed in the quark model 
to understand qa in neutron beta decay. 

Given these points, the extraction of jV^I from the measurements of T(B— >DeP e ), and 
T(B^D*eis e ) should be reliable. Ultimately, a precise determination of jV^I will come from 
a careful consideration of the heavy quark limit in which models like this one have been used 
to estimate l/m q and l/m Q corrections. 

The determination of \V u b\ is important for understanding CP violation in the standard 
model, since it is vital for determining the area of the"unitarity triangle" to which standard 
model CP violation is proportional. In most decays considered in this paper the model 
dependence of our results is modest. The b—>ueu e decays are, however, an exception. The 
large available recoil, the relativistic nature of the n and p, and the fact that such decays 
are far from any symmetry limits leave these predictions very exposed to uncertainties. We 
estimate that the theoretical uncertainties within our model for extracting from 
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the inclusive endpoint spectrum are at the 50% level. The uncertainties associated with 
individual exclusive channels are even larger: we would estimate them to be almost a factor 
of two for B — > ixtvt (where our model uncertainties are compounded by the uncertain 
effects of the nearby B* pole [46]) and 50% for B — > plve and B — > lu£u£. Fortunately, the 
determination of |T4&| can be greatly improved by combining the observation of exclusive 
B decays with their analogous D decays since such measurements can be related by Heavy 
Quark Symmetry. Here, once again, a model like ISGW2 has an important role to play, 
since it can assess the size of symmetry-breaking effects in this procedure. As we showed in 
Section V, our model predicts that this technique should allow the extraction of \V u b\ with 
a theoretical error of about 10%. 
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APPENDIX A: hyperfine-corrected wavefunctions 

Of the two leading order effects which break Heavy Quark Symmetry, the heavy quark 
kinetic energy and its hyperfine interaction, only the first was included in the constituent 
quark model which formed the basis of the ISGW prediction for form factors. In this 
Appendix we present a simple extension of their spectroscopic model which qualitatively 
takes hyperfine interactions into account. 

The extended spectroscopic model remains a nonrelativistic constituent quark model 
with essentially the same Coulomb-plus-linear central potential: 

V( r ) = + c + . (91) 

3r 

One cannot simply extend such a model by adding the non-relativistic hyperfine interaction 
since the Fermi spin-spin contact term, which is proportional to Si ■ Sj 5 3 (r) is an illegal 
operator in the Schrodinger equation. The problem is that, in channels where it is attrac- 
tive, this interaction is more singular than the kinetic energy so that the solutions of the 
Schrodinger equation collapse into = and to infinitely negative energies. This problem 
is solved by relativistic corrections which turn this operator into an extended and nonlocal 
one. Here we model this behavior by taking 

.ffijm,- i .32iTaa s Si ■ Sj 5 3 (f), r m i m j .i . . 

= < 92) 

where the term in parentheses would be the ordinary Fermi contact term if the anomalous 
coupling coefficient a were unity, and where Ei = (m 2 + p 2 )^. We have examined the effects 
of smearing out 5 3 (f) and found it to be small compared to the very strong nonlocality 
created by the pre- and postfactors of 5 . We have also ignored the tensor part of the 

hyperfine interaction as well as spin-orbit interactions. Neither can play a leading-order role 
in the S-waves which dominate our discussions, and both are also observed to be relatively 
weak even in excited states. Finally, we have made the running of a s in ISGW slightly less 
crude by taking a s for the p, K*, 0, D*, B, D*, B*, ip, B c , and T families to be 0.60, 0.55, 
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0.55, 0.50, 0.50, 0.45, 0.40, 0.40, 0.35, and 0.30, respectively on the basis of their reduced 
masses. Note that, following Ref. [35] we assume a s "freezes out" at 0.60 at low mass scales. 

Also as in ISGW, we solve the Hamiltonian variationally in a basis of harmonic oscillator 
states truncated to include only the IS, IP and 2S states. An interesting feature of this 
procedure, realized numerically in ISGW but thought to be an accident, is that the solution 
will exhibit zero IS — 2S mixing. This is proved in Appendix B. As a result, we need only 
the diagonal matrix elements of p 2 , 1/r, 1, r, and (Al) to solve for the wavefunctions. All 
but the last are trivial; for it one easily finds that the diagonal matrix elements all vanish 
except in the S- waves where 

(n 2s+1 S 2s+1 \H hyp \n^S 2s+1 ) = [ 2s{s + 1} ~ 3 ](g^)fe(0)| 2 (93) 

with s = or s = 1 the total quark-plus-antiquark spin, 

with mi and m 2 the constituent quark and antiquark masses and 4> n s{p) the momentum 
space wavefunction. (Note that ^„s(0) reduces to the nonrelativistic spatial wavefunction 
at f = in the nonrelativistic limit). 

On minimizing energies with respect to the gaussian wavefunction parameters $ previ- 
ously defined in ISGW and searching for a fit to the observed meson spectra, we found the 
results listed in Tables Al and A2. We assume that defects like the p — n splitting would 
improve with a larger basis space, but given the crude nature of this quasirelativistic model 
and our goal of a qualitative description of hyperfine effects, we do not attempt a better fit 
via a more complicated variant of the model. We also emphasize that full consistency would 
require a parallel relativistic treatment of both the spectrum and the weak matrix elements; 
this far more ambitious program would be very worthwhile, but it is well beyond the scope 
of this work. 
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APPENDIX B: variational solution in a 1S-2S basis 

Let [3^ be the value of P$ which minimizes the full Hamiltonian in the harmonic oscillator 



ground state 



Thus 



d/3i 



4f = (95) 



J d^r^f s s) H(p, r)# = (96) 



for p s = (3P. But 



= (h l/2 Ps l 4i s) (97) 



since 



dp s ri ° y 2< 



4r = (l) 1/2 ^-lpf)e-^ . (98) 



Thus, if pP minimizes H, 



J tnUpmlSP + / tntpHiUp = o . (99) 

But (V>is s) |#(p,r)lV4? S } ) = (4s s) \H(p,r)\i;[ p s s } }* and both are real so 

(tfrf) = o (ioo) 



i.e. ip2S does not mix with ipis if it has been chosen variationally to minimize H. This 
argument can be extended; for example, ip2P does not mix with ipip if it has been chosen 
variationally to minimize H. 
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APPENDIX C: form factor modifications for ISGW2 

As described in the text, the formulas of Appendix B of ISGW (and of those additional 
formulas in Refs. [25] and [30] needed when the lepton mass cannot be neglected) require 
modification. 

First, as discussed in Section III.C, all formulas are affected by the replacement in eq. 
(Bl) of ISGW shown in eq. (20). In addition, the conversion from the ff m to the described 
in Section III. A introduces factors of {fhs lm^ aB ^ a \fnx /fhx) nx ^ into each ISGW formula. 
With both of these changes effected, the factor F n of eq. (Bl) of ISGW is converted to a 
factor we denote by F^> since it now depends on the form factor a under consideration. The 
powers and %(a) required to make these conversions are given in Table CI. Note 

that in many instances it is a special combination {e.g., /+ + /_) which has a simple mass 
scaling law and not the individual form factors {e.g., f + ). For this reason we quote below 
formulas for these special combinations. To compute a particular form factor in such cases, 
one must apply the methods described here to those special combinations and then combine 
these results. Section III. A also describes how all S'-wave to S-wave transition form factors 
must be modified by the matching conditions given in Section III. A. These corrections lead 
below to the appearance of the factors = C^l + f3 {a) a s /n) { or, in the case of a + + a_, 
to terms proportional to Cjif3^a s /ir) which are the coefficients of £{w) in eq. (5). (In 
practice we use the "renormalization group improved" matching for all decays except those 
induced by the s — > u transition for which we implement "lowest order matching".) 

In the following we employ the notation of ISGW augmented with w as identified by Eq. 
(3). As in ISGW, all formulas are given for the b — > c transition but can be immediately 
adapted to any decay; see Ref. [1] for details including the explicit definition of each form 
factor. For an example, see eq. (103) below for the / form factor. The new formulas are: 

1. Eqs. (B8) and (B9) are replaced by the two equations 



/+ + /- 



2 - 



m x 



m q 



1 - 



m d m q f3 2 B 
2fi + m x f3 2 BX 



) 



F (f++f-) R (f++f-) 



(101) 
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/+-/- = 



m B 



\ m d m q , , 

v 2fi + m x p 2 BX J 3 



(102) 



which determine both / + and /_. 



2. Eq. (B15) becomes 

/ = C f m B 



[l + w) + 



m d (w - 1) 



(103) 



The (1 + w) and (w — 1) terms come from the constraints of Heavy Quark Symmetry 
in leading [4] and next-to- leading order [42-44], respectively, in the 1/toq expansion. 
Cf is from Eq. (7), and is the modified factor F 3 described above, namely 



<f) 



(mx/fhx) 2 (mB/m B ) 2 x (m x /m B ) 2 



i / PxPb 



P 2 



[l + ^r XB (t m -t)}- 



BX 



12 



(104) 



3. Eq. (B16) for g becomes 



1 



m d {3 B 



2fi_m x p BX 



F (9) R (g) 



(105) 



4. Eq. (B17) for a + is replaced by the two equations 



a + + a. 



m d 



m d (3 X 

[l + w)m q m b I3 2 BX V 2rh B (3 BX 



1 - 



7T 



m B 



i(a + +a_) 



(106) 



a + — a_ = — 



m d (3 2 x 
^+Plx 



wm d m B 



Px 



( m d P x 

(w + l)m q m b (3 BX \ 2m B (3 



1 - 



BX . 



p( a + ~ a - ) JjK a + ~ a - 



(107) 
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These formulas follow from ISGW and Ref. [25] with w-dependence dictated by the 
constraints of Heavy Quark Symmetry in order 1/mg [42-44]. They determine both 
a + and a_. 



5. Eq. (B23) for h becomes 



h 



m d 



1 



m d (3 B 



2fi_m x p BX 



(108) 



6. Eq. (B24) for k becomes 



k = 



m d 
V2f3 L 



[l + w)F£ 



(k) 



(109) 



7. Eq. (B25) for b + is replaced by the two equations 



b + + b_ 



ft 



i - 



m d (3 



A^2m q m h m B (3 B Pbx \ 2 ™b P 2 BX 



(110) 



6+ - 6_ = - 



m d 



V2m b m x f3 B 



1 - 



m d m b f3 2 x 
2fi + m B (3 2 BX 



| rn d (3 X 
4m q (3 2 BX 

which determine both b + and 6_. 



8. Eq. (B31) for q becomes 



m d 3 2 



x 



2m B (3 2 BX ^ 



<b+-b-) 



Q = ~ 



m d 



5 + w 
6 



2rh x p B 

See also eqs. (125) and (129) below. 

See ref. [23] of ref. [29] for an explanation of the sign change. 



9. Eq. (B32) becomes 
£ = -rh B /3 B 



1 m d m x {w — 1) / 5 + w 1 m d j3 



6m g 2/i_ mi 



(111) 



(112) 



(113) 



See also eqs. (122) and (126) below. 
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10. Eq. (B33) for c + is replaced by the two equations 



c + + c_ = 



m d m x 



2m q m B (3 B 



m d m q (3 
2m x /i„ (3 2 BX 



B \ t?(c++c-) 



(114) 



c + — c_ = — 



m d m x ( w + 2 m^m, /3 



^(c + -c_) 



2m g m B /3 B \ 3 2m x \i_ (3 BX 



(115) 



which determine both c + and c_. See also eqs. (123), (124), (127), and (128) below. 
11. Eq. (B37) for u+ is replaced by the two equations 



u + +U- — —\ 



'2rn d (u ++u _) 
3(3 B 5 



(116) 



u + — u_ 



l2 m d m B ^ u+ _ u _) 
3(3 B rh x 5 



(117) 



which determine both u + and u. 



12. Eq. (B43) for v becomes 



fh B (3 £ 



(w - 1) m d 
4:\/2mbm q fh x < o\[2 fn x (3 B 



+ 



(118) 



See also eqs. (125) and (129) below. 



13. Eq. (B44) for r becomes 



r = 



m B (3 £ 
V2 



1 m d m x _ 2 



(119) 



See also eqs. (122) and (126) below. 



14. Eq. (B45) for s + is replaced by the two equations 



m d 



V2m B /3 B V m q 2v+Pbx 



(120) 
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S + — S- = 



' 4-w _ m d m q (3% \ ( s+ s_) 
3 o™..,, ni I 5 



121) 



V2m q (3 B V 3 2mxfi+P 2 BX 
which determine both s + and s_. See also eqs. (123), (124), (127), and (128) below. 



15. Eq. (B49) for /' is replaced by the two equations 



/;+/: 



r^ L)u _rn ±v 



m, 



i J 



(122) 



/:-/' 




m x . 



,(/;-/!) 



(123) 



where 



Wlx 



+ 



3fe 



(124) 



and where 



V 



Q(3 2 



m 6 ' 



7- 



A 2 



(5 + r) 



(125) 



T =ml(3\{w-l) 



(126) 



These equations determine both /' and f_. 



16. Eq. (B55) for /' becomes 



f' = C f ^m B (l + w)UFP 



(127) 



where U is given above. 



17. Eq. (B16) for g' becomes 



9 



(— - 



m 



<9') 



(128) 



where once again U is given above. 
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18. Eq. (B17) for a' + is replaced by the two equations 

f <">';, j 



a' + a'_ = - 



PI 



3 m q m b p BX 



[1 + jr] 



3m 2 /% 



Sm B p B P BX 



+ 



3m d (3 x 



,(a' +o'_) 



(129) 



'2 3mg 
3 2mb7hx 



1 - 



[1 + ^] - 



m d # 2 



2m B/ 9 



1- 



5/% 
2 



3/3 



12m q m B (3 BX 



1- 



a 2 



+ 



BX 



ill 

7/3 2 



) *i 



(o' + -a'_) 



(130) 



BX 



These formulas determine both a' + and a'_. 

For transitions to excited P-wave heavy quark systems, heavy quark symmetry tells us 
that the L — S coupled states 3 Pi and 1 P 1 which are appropriate to the light 1=1 and 
1 = sectors and to states of definite C parity should be replaced by the j — j coupled 
states with sf = § and \ + . If we define form factors £3,c + 3, c_3, and q-i to be the exact 
analogs of £, c+, c_, and q of equations (B26) and (B27) of ISGW but for the sj e = | + state 
with J p = 1 + , and a parallel set £1, c, i,c_i, and 01 for the sl l = \ + state with J p = 1 + , 

2 " r 2 2 2 * 2 

then 



£3 = 
2 



,0 



V3 



+ 



m x mAw 



w 



m d f3_ 



2m q 2/j,_m x p BX 



(<a) 



(131) 



V3m d 

C 1 3 + C_3 = — — - — - — 

5 2p B m B 



1 - 



3m„ 



m d p 2 
3(3 2 



B 



BX 



'—--) 



(c, 3 +C 3 ) 
F + 3? -1 



(132) 
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C i 3 — C_3 
^2 2 



2V3/3. 



(2 - w)m x m d /? s 



P. 



BX 



2/z_ 



(c,3-c 3) 
F + 2 -7 



(133) 



and 



2^3 \ 2 2m d m q m h J ftm^ 5 



(134) 



2m„ 



2m ; 



+ 



m d fh x {w - 1) 



m 9 2fi_m x (3 B x_ 



(ti) 



(135) 



c.i 
^2 



2 



™2 /32 



VQrn B m q (3 B (3 BX 



( c +i+ c -i) 



(136) 



c , 1 — c 

"•"2 



/2 m d 
3rhx(3 B 



1 + 



2m q (3 2 BX 



(137) 



9i = 



1 



(138) 



. 6 { Am d m q m b J f3 B m x 
We use these latter formulas for decays to all heavy-light final states, including the kaons. 



Of course, to the extent that the s 



3 + 



and multiplets are degenerate there formulas 



will give total rates to the two 1 + states that are identical to the 3 Pi and l Pi formulas. 
However, the latter are needed to predict the rates to individual states, e.g., the rate for 

b -> r>r W 
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Table I. Relativistic corrections to the / form factor 



decay 


C f for IS 


C f for 2S 


D — > p, LO 


0.889 


0.740 


D -> K* 


0.928 


0.782 


D S ^K* 


0.873 


0.739 


D s ^4> 


0.911 


0.773 


B — ► p, to 


0.905 


0.776 


B ^D* 


0.989 


0.929 


B S ^K* 


0.892 


0.781 


B s - D* s 


0.984 


0.924 


B C ^D* 


0.868 


0.779 


B c ^^ 


0.967 


0.899 
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Table II. Exclusive partial widths for the b^c semileptonic decays, B— >X c jeu e , B s ^X cS ei? e 
and B c — > X cg ez/ e , in units of 10 13 | Vb c \ 2 sec -1 . The Heavy Quark Symmetry notation rf l Lj 
is used for the final states with unequal mass quarks b \ Also included are the physical 
meson masses used (in GeV), taken from Ref. [33] if possible; properties of unobserved 
or controversial states (given in parentheses) are taken from Ref. [35]. The masses of the 
decaying particles (in GeV) are 5.28, 5.38 d \ and (6.27). 







Bd^X cd ev e 




B s ^X C geu e 


B c 


— »■ X c5 eu e 


X 


mass 


partial width 


mass 


partial width 


mass 


partial width 


1 \ <?n 


1.87 


1.42 


1.97 


1.31 


2.98 


0.99 


l*Si 


2.01 


2.81 


2.11 


2.49 


3.10 


1.57 




2.46 


0.10 


2.57 c ) 


0.14 


3.56 


0.12 


liPi 


2.42 


0.20 


(2.54) 


0.25 


3.52 fe ) 


0.21 


l3Pl 


(2.49) 


0.04 


(2.57) 


0.05 


3.51 fe ) 


0.08 


l^Po 


(2.40) 


0.03 


(2.48) 


0.04 


3.42 


0.04 




(2.58) 


0.00 


(2.67) 


0.01 


(3.62) 


0.07 


22^1 


(2.64) 


0.06 


(2.73) 


0.13 


3.69 


0.29 


total 


4.66 


4.41 


3.36 



a ) The Tl/Tt values for these decays are 1.08, 1.03, and 0.88, respectively. 

b ) We list B c -> ^Pi (i.e., xa) under l§Pi and B c -> l 3 Pi (i.e., h cl ) under \\p x . 

c ) See Ref. [36]. 

d ) See Ref. [37]. 
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Table III. Exclusive partial widths for the c^s semileptonic decays, D— >X S ue + v e , 
D s ^X sS e + u e and B c ^X s ie + u ei in units of 10 10 | V cs \ 2 sec -1 . The Heavy Quark Symmetry 
notation n Se Lj is used for the final states with unequal mass quarks c \ Also included are 
the physical meson masses used (in GeV), taken from Ref.[33] if possible; properties of un- 
observed or controversial states (given in parentheses) are taken from Ref . [35] . The masses 
of the decaying particles (in GeV) are 1.87, 1.97, and (6.27), respectively. 





D 




D s 


~\s~ T- 


B c 




X 


mass 


partial width 


mass 


partial width 


mass 


partial width 


1 


0.49 


10.5 


0.55 6 ) 
0.96 6 ) 


3.7 
3.2 


5.38 d ) 


2.2 


US^ 


0.89 


5.7 


1.02 


4.8 


(5.45) 


2.7 


lip 2 


1.43 


0.00 


1.53 


0.00 


(5.88) 


0.00 


l§Pi 


1.27 


0.34 


1.38 c ) 


0.27 


(5.88) 


0.06 


l5Pl 


1.40 


0.00 


1.51 c ) 


0.03 


(5.88) 


0.00 


15P 


1.43 


0.00 


1.52 


0.00 


(5.88) 


0.00 


25 S 


(1.45) 


0.00 


(1.63) 


0.00 


(5.98) 


0.01 


25 Si 


(1.58) 


0.00 


(1.69) 


0.01 


(6.01) 


0.01 


total 


16.6 


12.1 


5.0 



a ) The Yl/Yt values for these decays are 0.94, 0.96, and 1.03, respectively. 

fe ) We use the approximation of ideal mixing in I — states in every sector except the 
ground state pseudoscalars where we assume an r\ — r( mixing angle of —20°. If this mixing 
angle were changed to —10°, then the entries in the Table to rj and rj' would change to 5.6 
and 2.4; note that while the individual rates change substantially, the total rate to these 
two states would only increase by about 16%. 

c ) For D s — > X sS we list the rate to the l x Pi under laP 1 and that for the l 3 Pi state under 
lsPi. 

^ See Ref. [37]. 
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Table IV. Exclusive partial widths for the c—*d semileptonic decays D — >X dii e + v e and 
D s ^X ds e + v e) in units of 10 10 1 X^rf| 2 sec -1 . Also included are the physical meson masses used 
(in GeV), taken from Ref. [33] if possible; properties of unobserved or controversial states 
(given in parentheses) are taken from Ref. [35]. The masses of the decaying particles (in 
GeV) are 1.87 and 1.97, respectively. 







D°^X m e + u e 






X 


mass 


partial width 


mass 


partial width 




0.14 


9.8 


0.50 


8.9 


l3c a) 
i o\ 


0.77 


4.9 


0.89 


4.4 


i 3 p 2 


1.32 


0.01 


1.43 


0.01 


i 1 ^ 


1.23 


0.52 


1.27 6 ) 


1.5 


i 3 Pi 


1.26 


0.32 


1.40 6 ) 


0.01 


i 3 p 


1.30 


0.00 


1.43 


0.00 




1.30 


0.02 


(1.45) 


0.04 




(1.45) 


0.03 


(1.58) 


0.03 


total 


15.6 


14.9 



a ) The T L /T T values for these decays are 0.67 and 0.76, respectively. 

fe )We list the rate to the mainly l^Pi state under I 1 Pi and that for the mainly lap 1 under 
l 3 Pi state. 
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Table V. Exclusive partial widths for the c-^d semileptonic decay D + ^X dd e + v e in units of 
10 10 1 V cd \ 2 sec -1 , separated into 1 = 1 and 1 = contributions. Also included are the physical 
meson masses used (in GeV), taken from Ref. [33] if possible; properties of unobserved or 
controversial states (given in parentheses) are taken from Ref. [35] . The mass of the decaying 
particle is 1.87 GeV. 







D + ^X dd -e + v e 




D + ^X dd -e + v e 






1=1 




1 = 


X 


mass 


partial width 


mass 


partial width 


1 1 a 


0.14 


4.9 


0.55 fc ) 


3.0 








0.96 6 ) 


0.6 




0.77 


2.5 


0.78 


2.4 


1 3 P 2 


1.32 


0.00 


1.27 


0.00 




1.23 


0.26 


1.17 


0.39 


1 3 P1 


1.26 


0.16 


1.28 


0.13 


1 3 P 


1.30 


0.00 


1.30 


0.00 


2% 


1.30 


0.01 


(1.44) 


0.00 


2 3 5i 


(1.45) 


0.02 


(1.46) 


0.01 


total 


7.8 


6.5 



a ) The Tl/Yt values for these decays are 0.67 and 0.68, respectively. 

fe ) We use the approximation of ideal mixing in I — states in every sector except the 
ground state pseudoscalars where we assume an 77 — 77' mixing angle of —20°. If this mixing 
angle were changed to —10°, then the entries in the Table to i] and rj' would change to 2.2 
and 0.9; note that while the individual rates change substantially, the total rate to these 
two states would only decrease by about 14%. 
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Table VI. Exclusive partial widths for the b— m semileptonic decays with a light spectator, 
£>°^X u jez/ e , and B~ ^X uU eu e , in units of 10 13 | V^| 2 sec" 1 . Also included are the physical 
meson masses used (in GeV), taken from Ref. [33] if possible; properties of unobserved or 
controversial states (given in parentheses) are taken from Ref. [35]. (The masses for the 
1 = 1 final states in B~^X ua eis e are the same as those for B Q -^X u ^ev e ) The mass of the 
decaying particle is 5.28 GeV. 











1=1 


X uu &v e 


1 = 


X 


mass 


partial width 


mass 


partial width 


mass 


partial width 




0.14 


0.96 


0.14 


0.48 


0.55 6 ) 
0.96 6 ) 


0.45 
0.28 


1 3 S1 a ) 


0.77 


1.42 


0.77 


0.71 


0.78 


0.71 


1 3 P 2 


1.32 


0.33 


1.32 


0.16 


1.27 


0.18 


I 1 Pi 


1.23 


1.09 


1.23 


0.54 


1.17 


0.57 


1 3 P1 


1.26 


0.87 


1.26 


0.43 


1.28 


0.41 


l 3 Po 


1.30 


0.05 


1.30 


0.02 


1.30 


0.03 


2 1 S 


1.30 


0.17 


1.30 


0.08 


(1.44) 


0.08 


2 3 5i 


(1.45) 


0.41 


(1.45) 


0.20 


(1.46) 


0.20 


partial total 


5.3 




2.6 




2.9 



a ) The T L /T T values for these decays are 0.30, 0.30, and 0.30 respectively. 

^ We use the approximation of ideal mixing in I — states in every sector except the 
ground state pseudoscalars where we assume an rj — r( mixing angle of —20°. If this mixing 
angle were changed to —10°, then the entries in the Table to rj and rj' would change to 0.34 
and 0.41; note that while the individual rates change substantially, the total rate to these 
two states would only increase by about 3%. 
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Table VII. Exclusive partial widths for the b— m semileptonic decays with a heavy spectator, 
B s ^X uS eu ei and B c — > X u5 eu e , in units of 10 13 | sec -1 . The Heavy Quark Symmetry 
notation n Se Lj is used for the final states with unequal mass quarks. Also included are 
the physical meson masses used (in GeV), taken from Ref. [33] if possible; properties of 
unobserved or controversial states (given in parentheses) are taken from Ref. [35]. The 
masses of the decaying particles (in GeV) are 5.38 6 ^ and (6.27), respectively. 







B s ^X uS ev e 




B c — > X U cev e 


X 


mass 


partial width 


mass 


partial width 


1 * >->o 


0.49 


0.85 


1.87 


0.30 


13^1 a ) 


0.89 


1.14 


2.01 


0.62 


lip 2 


1.43 


0.28 


2.46 


0.06 


lfPi 


1.27 


1.72 


2.42 


0.62 


l3Pl 


1.40 


0.08 


(2.49) 


0.04 


HPo 


1.43 


0.04 


(2.40) 


0.01 




(1.45) 


0.45 


(2.58) 


0.46 


25^1 


(1.58) 


0.54 


(2.64) 


0.40 


partial total 


5.1 


2.5 



a ) The Yl/Y t values for these decays are 0.45 and 0.61, respectively. 
6 WRef. [37] 
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Table VIII. Comparison at zero recoil of the ISGW2 meson form factors without pertur- 
bative matching corrections to those of leading order Heavy Quark Symmetry and Heavy 
Quark Symmetry including the Oil/mq) corrections predicted by ISGW2. Note that these 
ISGW2 form factors cannot be directly compared to experiment: they are just the /J m scaled 
by factors of (rh B /m B ) nB ^ (rh x /rn x ) nx ^ with the exponents of Table CI. The form factors 
are those for mesons corresponding to a light (u or d) spectator. 





h 


f- 


9 


/ 


(a+ + a_) 


(a+ — a_) 


ISGW2-no matching 


b^cev e 


1.00 


-0.03 


1.11 


0.97 


-0.08 


1.04 


ISGW2-no matching 




0.98 


-0.00 


1.28 


0.84 


-0.23 


1.18 


HQS 




1 





1 


1 





1 


HQS 




1 





1 


1 





1 


HQS+o(^) 


b^>cei> e 


1 


-0.06 


1.12 


1 


-0.09 


1.03 


HQS+o(^) 




1 


-0.21 


1.39 


1 


-0.30 


1.09 
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Table IX. Predictions for the six form factors for B -> Dlv t and B -> D*^. The HQET 
column shows the effect of QCD radiative corrections alone to the HQS symmetry limit 
column. 





TSOW2 

1 kJ V.I VV £d 


mow 

1 >J V I VV 


HOS 




f-^ (t>rn) 


1.00 


1.01 


1 


1.00 


f-(tm) 


-0.09 


-0.05 





-0.02 


g{t m ) 


1.17 


1.12 


1 


1.06 




0.91 


1.00 


1 


0.93 


2a + (t m ) 


0.83 


0.95 


1 


0.88 


2a_ (t m ) 


-1.19 


-1.11 


-1 


-1.08 
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Table X. Comparison of the form factors for D-^Ke + v e and D^>K*e + is e with experiment. 
We have used the t dependence assumed in the fits to data to extrapolate the experimental 
form factors to t — t m from t — 0. 





experiment [23] 


ISGW2 


ISGW 


f+{tm) 


1.42 ±0.25 


1.23 


1.16 


f(t m ) (GeV) 


2.21 ±0.19 


1.92 


2.76 


g(t m ) (GeV- 1 ) 


0.55 ± 0.08 


0.55 


0.47 


a+(t m ) (GeV- 1 ) 


-0.21 ±0.04 


-0.34 


-0.37 
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Table Al. quark model parameters 



parameter 






b 


0.18 GeV^ 


0.18 GeV 2 


c 


-0.81 GeV 


-0.84 GeV 


a s 


0.60 ->■ 0.30 (see text) 


0.50 ->■ 0.30 (see ISGW) 


m u = m d 


0.33 GeV 


0.33 GeV 


1 1 b S 


55 GeV 


55 GeV 


m c 


1.82 GeV 


1.82 GeV 


m b 


5.20 GeV 


5.12 GeV 


a 


2.8 


not applicable 
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Table A2. the masses and (3 values in GeV for variational solutions of the hyperfine- 
corrected Coulomb plus linear problem in the IS", IP, 2S basis 





meson flavor: 


ud 


us 


ss 


cu 


cs 


ub 


sb 


cc 


be 




m 


0.35 


0.55 


0.62 


1.86 


1.94 


5.27 


5.33 


2.95 


6.33 


VS 



























0.41 


0.44 


0.53 


0.45 


0.56 


0.43 


0.54 


0.88 


0.92 




m 


0.74 


0.87 


0.97 


2.01 


2.10 


5.33 


5.40 


3.13 


6.42 


l 6 Si 
























P 


0.30 


0.33 


0.37 


0.38 


0.44 


0.40 


0.49 


0.62 


0.75 




m 


1.24 


1.35 


1.42 


2.48 


2.53 


5.81 


5.84 


3.51 


6.79 




























0.28 


0.30 


0.33 


0.33 


0.38 


0.35 


0.41 


0.52 


0.60 
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Table CI: The factors F n of eq. (Bl) of ISGW are to be replaced by F^' which have 
the modification shown in eq. (27) of Section III.C and are multiplied by the a-dependent 

/ _ snj(a) / - \i»i(a) 

factor with n^a) and nx{ot) given here. 



TABLE I. 



form factors 






(a) 


n B (a) 


n x (a) 


/+ + /-,/; + /: 


-1/2 


+1/2 


f+ — f-> f+ — f- 


+1/2 


-1/2 


9,9' 


-1/2 


-1/2 


/,/' 


+1/2 


+1/2 


a + + a_ , a' + + a'_ 


-3/2 


+1/2 


a + — a_ , a' + — a'_ 


-1/2 


-1/2 


/i 


-3/2 


-1/2 


fc 


-1/2 


+1/2 




-5/2 


+1/2 


JL T 

o + — o_ 


—3/2 


1 /a 

-1/2 


2 2 


+1/2 


+1/2 


C+ + C_ , S + + S_ , C , 3 + C 3 , C , 1 + C 1 

^2 2^2 2 


-3/2 


+1/2 


C+ — C- , S+ — S- , C , 3 — C 3 , c , 1 — C 1 

" r 2 2^2 2 


-1/2 


-1/2 


q,v,qs,qi 
2 2 


-1/2 


-1/2 


u + + u_ 


-1/2 


+1/2 


u + — u_ 


+1/2 


-1/2 
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Fig. 1. Measurements of the pion form factor [18] compared to the form factors of 
ISGW and ISGW2. 
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Fig. 2. (1/T)( dT / dE e ) for i? — >■ X c ez/ e showing contributions of D } D* , and the total 
contribution from all IS, IP, and 2S states; also shown is the corresponding free quark curve. 
Absolute rates are given in Table II, and may be compared to rj ree = 4.84 X 1 13 1 V^j, | 2 sec -1 . 
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Fig. 3. (1/T)( dT / dE e ) for B s — > X C gev e showing contributions of D S} D s * } and the total 
contribution from all IS, IP, and 2S state s; also shown is the corresponding free quark curve. 
Absolute rates are given in Table II, and may be compared to T j ree = 4.84 X 1 13 1 V^j, | 2 sec -1 . 
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Fig. 4. (l/T)(dT / dE e ) for B c — > X C cev e showing contributions of rj c , J/^, and the total 
contribution from all IS, IP, and 2S states; also shown is the corresponding free quark curve. 
Absolute rates are given in Table II, and may be compared to rj ree = 4.84 X 1 13 1 V^j, | 2 sec -1 . 
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Fig. 5. (1/T)( dT / dE e ) for D — > X sU e + p e showing contributions of K, K* } and the total 
contribution from all IS, IP, and 2S states; also shown is the corresponding free quark curve. 
Absolute rates are given in Table III, and may be compared to rj ree = 0.34 X 1 12 1 V^ s | 2 sec -1 . 
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Fig. 6. (1/T)( dT / dE e ) for D — > X S ge + v e showing contributions of rj } rj' } J/^, and 
the total contribution from all IS, IP, and 2S states; also shown is the corresponding free 
quark curve. Absolute rates are given in Table III, and may be compared to rj ree = 
0.34 x lO^lKspsec" 1 . 
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Fig. 7. (l/T)(dT / dE e ) for B — > X s ie + v e showing contributions of B S} B s * } and the total 
contribu tion from all IS, IP, and 2S states; also shown is the corresponding free quark curve. 
Absolute rates are given in Table III, and may be compared to rj ree = 0.34 X 1 12 1 V^ s | 2 sec -1 . 
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Fig. 8a. (l/T)(dT / dE e ) for D° — > X^e + z/ e showing contributions of 7r, p } and the total 
contribution from all IS, IP, and 2S states. Absolute rates are given in Table IV. 
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Fig. 8b. (l/T)(dT / dE e ) for D + — > X dd e + v e showing contributions from 7r°, rj } ?/, p and u, 
and the total including IP and 2S states. Absolute rates are given in Table V. The free quark 
decay curve is shown on both plots and has an absolute rate of rf ree = 0.51 X 10 12 1 | 2 sec -1 
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Fig. 9. (l/T)(dT / dE e ) for D s — s- X<ise + v e showing contributions of K° } K*°, and the total 
contribution from all IS, IP, and 2S states; also shown is the corresponding free quark curve. 
Absolute rates are given in Table IV, and may be compared to rf ree = 0.51 X 10 12 \V c d\ 2 sec -1 . 
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Fig. 10a. (l/rf ree )( dT / dE e ) for B° — > X u ^ev e showing contributions of 7r, p } the 
states, and the 2S states tt' and p'. 
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Fig. 10b. (l/rf ree )( dT / dE e ) for B~ — > X uU ev e showing contributions of 7r°, 77, ?/, p° } lo } 
and the IP and 2S states. The corresponding free quark curve is shown for both graphs 
corresponding to rj ree = 1.15 X 10 14 |Vf, u | 2 sec -1 . The partial widths are given in Table VI. 
Note that the curves shown are all normalized to rf ree since our partial sum over exclusive 
channels does not exhaust the semileptonic rate. 
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Fig. 11. (l/rf ree )( dT / dE e ) for B s — > X U gev e showing contributions of K } K* } the IP 
states, and the 2S states; also shown is the corresponding free quark curve. Absolute rates 
are given in Table VII, and may be compared to rf ree = 1.15 X 10 14 |V6„| 2 sec -1 . 
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Fig. 12. (l/rf ree )( dT / dE e ) for B c — > X U cev e showing contributions of D } D* , the IP 
states, and the 2S states; also shown is the corresponding free quark curve. Absolute rates 
are given in Table VII, and may be compared to rf ree = 1.15 X 10 14 |V6„| 2 sec -1 . 
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